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Abstract 

The universe is smooth on large scales but very inhomogeneous on small scales. Why is the 
spacetime on large scales modeled to a good approximation by the Friedmann equations? Are we 
sure that small-scale non-linearities do not induce a large backreaction? Related to this, what is 
the effective theory that describes the universe on large scales? In this paper we make progress 
in addressing these questions. We show that the effective theory for the long-wavelength universe 
behaves as a viscous fluid coupled to gravity: integrating out short-wavelength perturbations 
renormalizes the homogeneous background and introduces dissipative dynamics into the evolution 
of long- wavelength perturbations. The effective fluid has small perturbations and is characterized 
by a few parameters like an equation of state, a sound speed and a viscosity parameter. These 
parameters can be matched to numerical simulations or fitted from observations. We find that the 
backreaction of small-scale non-linearities is very small, being suppressed by the large hierarchy 
between the scale of non-linearities and the horizon scale. The effective pressure of the fluid is 
always positive and much too small to significantly affect the background evolution. Moreover, 
we prove that virialized scales decouple completely from the large-scale dynamics, at all orders 
in the post-Newtonian expansion. We propose that our effective theory be used to formulate 
a well-defined and controlled alternative to conventional perturbation theory, and we discuss 
possible observational applications. Finally, our way of reformulating results in second-order 
perturbation theory in terms of a long-wavelength effective fluid provides the opportunity to 
understand non-linear effects in a simple and physically intuitive way. 
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1 Introduction and Summary 



The formation of large-scale structures in the universe is characterized by two fundamental scales: 

i) the (comoving) Hubble scale 1-L~ l {ri) defines the extent of the observable universe at any 
given time r\ and limits the range over which interactions can influence the evolution of 
large-scale perturbations; 

ii) the non-linear scale k^(j]) describes the size of structures whose density contrast <5 (77, a;) = 
[p(rj,x)/p(rj) — 1] exceeds unity. 

The large hierarchy between these two scales, e = k^/W' 1 <C 1, is of fundamental importance for 
cosmology. It is this hierarchy of scales that is responsible for the success of linear perturbation 
theory: the most important features of the anisotropics of the cosmic microwave background 
(CMB) and the large-scale structure (LSS) observed in galaxy surveys are accurately described 
by linear perturbations around a homogeneous Friedmann- Robertson- Walker (FRW) background. 
However, with the advance of observations the study of small non-linear corrections to the long- 
wavelength dynamics is becoming more and more relevant. 

UV-IR coupling in cosmology. While for linearized perturbations different Fourier modes evolve 
independently, at the non-linear level two short-wavelength (UV) perturbations can couple to 
produce a long- wavelength (IR) perturbation, i.e. starting at quadratic order Fourier modes 
don't evolve independently. Beyond linear perturbation theory short-wavelength perturbations 
can therefore, in principle, affect the evolution of the long-distance universe. In particular, small- 
scale, non-linear fluctuations can give subtle backreaction effects both on the evolution of the 
background spacetime and the dynamics of long-wavelength perturbations. In fact, it has been 
claimed that the renormalization of the background from small-scale structures can be large 
enough to explain the acceleration of the universe without the need for dark energy [1]. In 
addition, the effect of small-scale non-linearities on the evolution of superhorizon perturbations 
from inflation provides an interesting case study [2] . These effects and others will be explored in 
this paper. 

We expect corrections to the linear evolution of scales with wavelengths comparable to (or 
larger than) the Hubble scale to be suppressed by the aforementioned hierarchy between the 
non-linear scale and the Hubble scale. This decoupling of short-wavelength (high-energy) fluctu- 
ations from the long-wavelength (low-energy) theory is of course a common feature of effective 
field theories [3] (see e.g. [4, 5, 6, 7] for recent reviews). It should therefore be possible to derive 
a long-wavelength effective theory in which corrections to the linear evolution appear system- 
atically suppressed by powers of e 2 (corrections linear in e are forbidden by the isotropy of the 
background). In this paper we aim to formalize the effective field theory approach as applied 
to cosmological perturbations. Specifically, we are interested in an effective description of the 

'in a universe filled with a combination of different matter components there are additional scales associated 
with the horizon at the time of transition between the different eras; e.g. in our universe the horizon at matter- 
radiation equality defines the scale = 'H _1 (?7eq)- This complication is of little consequence for the arguments 
made in this section, but will of course be taken into account in the remainder of the paper. 
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long-wavelength universe obtained by 'integrating out' short-wavelength modes. 2 We will derive 
the small corrections to the effective theory of long-wavelength perturbations arising from the 
UV-IR coupling imposed by the non-linearities of the Einstein equations and the matter sources. 

Matter fluctuations and perturbation theory. When following this logic, one may worry that the 
density contrast 5 becomes large below the non-linear scale and that small scales therefore lead 
to large backreaction effects on the long-wavelength modes. However, while the density contrast 
indeed becomes large, the spacetime perturbations and the particle velocities remain small (at 
least outside of black holes). The system is therefore still amenable to perturbation theory if 
organized in terms of the gravitational potential and the average particle velocity v rather 
than the density perturbation 5. Such an analysis reveals that very small scales in fact decouple 
from the large-scale evolution. We believe that this decoupling of short-wavelength non-linearities 
should even apply in the extreme case that the universe is filled with a gas of black holes. In this 
case, our perturbative scheme breaks down, but the effective theory can be matched continuously 
to the effective theory for the dynamics of black holes by Goldberger and Rothstein [8] , making 
a large backreaction of gravitational non-linearities even in this case unlikely. 

Below the virial scale there exist definitive relations between the gravitational potentials 
<I> and the velocities v, with large cancellations between the potential and kinetic energies. One of 
the main results of our work will be a proof that virial scales indeed decouple completely - i.e. they 
don't even lead to e 2 suppressed contributions to the effective pressure, although they of course 
lead to a small renormalization of the background density. In other words, the non-linear source 
terms for the evolution of large-scale modes can be expressed in a form that vanishes identically 
in the virial limit. We stress that this is more than standard effective field theory decoupling: 
indeed, according to the latter, the leading long-distance effect of short-distance physics is just 
a renormalization of the parameters of the long-distance effective theory [3]. Here instead, we 
claim that one such parameter - the effective pressure - does not even get renormalized in the 
virial limit. Therefore, our result goes more properly under the name of a 'non-renormalization 
theorem'. This significantly constrains our expectations for backreaction effects of small-scale 
non-linearities on the long-distance cosmological dynamics. We stress that the decoupling of 
virialized structures holds at all orders in the post-Newtonian expansion. Therefore it applies 
equally well to relativistic virialized systems, like for instance those containing black holes. 

The effective theory. Ultimately, our theory will be formulated as an FRW universe with 
small (quasi-linear) long-wavelength perturbations evolving in the presence of an effective fluid 
whose properties are determined by non-linear short-wavelength modes. We present the details 
of the effective theory in §5. The key element of the theory is the effective stress-energy tensor 
Tp V induced by the short-wavelength modes. Given its importance 3 , we will derive t^ v in two 

2 We define this procedure in detail in §4. There we explain that integrating out short-wavelength fluctuations 
amounts to smoothing the equations of motion and taking expectation values of the short-wavelength modes in the 
presence of long-wavelength perturbations, so that one is left with equations in terms only of the long-wavelength 
modes. 

3 The non-linear scalar source terms captured by r M „ have the following effects: 

1. the generation of vector perturbations o; t [9, 10, 11]; 

2. the generation of tensor perturbations Xif [H> 12, 13, 14]; 
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different ways: 

Effective stress-energy via (post-)Newton - constructive approach. We aim at understanding 
the effects of the short-scale non-linearities on the background Hubble expansion and on the 
long-wavelength fluctuations. In §3 we will study non-linear cosmological perturbations in a 
general-relativistic framework. However, we will encounter a number of technical complications 
that hide to some extent the physical relevance and intuitive nature of our findings. Fortunately, 
there are two important points emerging from the former discussion that simplify the problem 
considerably, allowing us to give a quicker yet rigorous derivation. First, the scale at which 
non-linearities in the perturbation equations become relevant is much smaller than the horizon 
scale. This allows us to concentrate on subhorizon distances and neglect the general-relativistic 
effects associated with the background FRW expansion. Second, for non-relativistic structures 
like clusters or galaxies, these non-linearities involve the matter sector only. That is, the short- 
scale gravitational dynamics is Newtonian to a very good approximation. Of course, the effect 
that these non-linear structures then have on the long-scale perturbations is posi-Newtonian in 
nature - since it involves the coupling of gravity to itself - but given the above considerations 
the following simplified approach suggests itself: short-scale perturbations evolve according to 
flat-space Newtonian equations, where all gravitational fields (short-scale as well as long-scale) 
are encoded in the Newtonian potential <£. The total stress-energy tensor t^ u of this system is 
conserved, in the ordinary sense, i.e. d^T^ = 0. Indeed, the existence of a conserved stress- 
energy tensor follows from locality and from invariance under space-time translations, regardless 
of Lorentz-invariance. The tensor t^ v has a gravitational contribution, of order p&, because the 
gravitational potential energy obviously participates in the stress-energy conservation already 
in Newtonian physics. We then smooth this stress-energy tensor over some scale A -1 larger 
than the typical inhomogeneity scale, and define an effective long-scale stress-energy tensor. The 
post-Newtonian leap is now to declare that the effective stress-energy tensor thus defined is 
what couples to long-wavelength gravitational fields. However, there is not much freedom in this 
assumption. Long- wavelength gravitational fields must couple to anything (including themselves) 
through a locally conserved symmetric tensor [15, 16]. For any given system, the only tensor with 
such properties is the (symmetrized) stress-energy tensor, which is unique up to total derivative 
terms of the form [16] 

dadpYpdW . (1.1) 

Here, the tensor £ is symmetric under the exchange of the two index pairs, and antisymmet- 
ric within each pair. The addition of a term of the form of (1.1) to a system's stress-energy 
tensor is identically conserved, does not affect the associated global charges (i.e., the total four- 
momentum), and most importantly for our purposes vanishes like two powers of momentum in 
the small- momentum limit. This ambiguity thus belongs in the class of 'higher-derivative' cor- 
rections (which we will discuss below), and can be neglected for long wavelengths. We work out 
the details of this approach in §2. 

3. the superhorizon evolution of scalar perturbations & 2 - ) 7^ [2] (see §6.4); 

4. the viscous damping of density perturbations (see §6.2). 
All of these effects are absent in linear perturbation theory. 
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Effective stress-energy via Einstein - deductive approach. A second equivalent description of the 
UV-IR coupling of cosmological fluctuations arises from a simple reorganization of the Einstein 
equations. First, we decompose the Einstein tensor into a homogeneous background (denoted 
by overbars) and terms that are linear (L) and non-linear (NL) in the metric perturbations, 
collectively denoted by SX(rj,x) = X(n,x) — X{rj). The Einstein equations can then be written 
as 

G^[X] + (G^) L [5X] + (G^) NL [SX 2 } = 8rfT p . (1.2) 

The background equations, = 8-kGT^, and the linearized Einstein equations, (G^) 1 = 
8irG (T M j,) L , are then defined in the standard way. The non-linear Einstein equations can be 
written in a form that is very similar to the linear equations, 

{G^) l = 8ttG{t^-T^), (1.3) 

where we defined the effective stress-energy pseudo-tensor 

r M , = ^- ( %^ L . (1.4) 



The stress-energy pseudo-tensor has a long history in studies of General Relativity (e.g. as an 
approach to studying gravitational raditation; see the books by Weinberg [18] or Landau and 
Lifshitz [19]) and will play a key role in this paper. 

Properties of the effective fluid. Given the form of t^ v in terms of short-scale (high-momentum) 
fields - <& s , v s - we can analyze its effects on the large-scale modes - 5g, $£, V£ - and on the 
homogeneous background. For the benefit of the impatient (curious) reader, we will now state 
some of the highlights of that analysis, leaving detailed derivations and explanations to the main 
text and the appendices: 

1. In the absence of long- wavelength perturbations or on very large (superhorizon) scales, 
the gravitational small-scale (subhorizon) non-linearities mimic an isotropic fluid whose 
effective density and pressure simply renormalize the background by terms of order of the 
velocity dispersion, (vg). The effective pressure of the fluid is always positive and much 
too small to significantly affect the background evolution. Moreover, the spatial part of the 
stress-energy tensor is equal to the second time-derivative of the moment of inertia tensor 

where [. . . ]a denotes spatial averaging over a region of size A -1 , and Jy is the moment 
of inertia associated with the same region. This shows that virialized structures decouple 
completely from the effective theory at large scales. The backreaction effects that we 
capture in our effective treatment therefore cannot explain the acceleration of the universe. 
Finally, the small induced pressure and the associated renormalization of the background 
explain the apparent superhorizon evolution of primordial curvature perturbation ( [2]. 
After the renormalization of the background is taken into account, £ is indeed constant on 
superhorizon scales. 
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2. The fluid is an imperfect fluid in the sense that the small-scale non-linearities induce dis- 
sipative terms and non-negligible anisotropic stress into the evolution of long-wavelength 
perturbations. At long wavelengths the fluid is characterized by only a few parameters like 
an equation of state, a sound speed and a viscosity parameter. 4 For instance, to leading 
order, the source term in the Euler equation may be written as 

hkj ([Tjj a) _ ji x 2 kjv\ 

-^^—-^-c vis —, (1.6) 

where c s and c v ; s are time-dependent coefficients. These parameters can be calibrated by 
computing the small-scale dynamics with numerical iV-body simulations. Alternatively, 
the fluid parameters may simply be retained as free parameters to be measured by fitting 
predictions of the effective theory to observations. This kind of matching calculation is of 
course common in effective field theory. 

3. The short- wavelength fluctuations provide a source of noise to the dynamics of the long- 
wavelength perturbations. Although this statistical noise has a negligible effect on the 
evolution of the background cosmology, it is not irrelevant in all contexts. For example, 
the stochastic contribution to the pressure fluctuations can be comparable to the pressure 
fluctuations induced by long-wavelength density fluctuations 6g for a wide range of scales. 
However, these fluctuations are uncorrelated at leading order with the amplitude of the 
long-wavelength modes, which results in a suppression of their importance in averaged 
quantities such as the power spectrum. 

An alternative to conventional perturbation theory. After the fluid parameters are determined 
from iV-body simulations of scales with high momenta, k > A, the effective fluid has small 
expansion parameters - {5e,$e,ve} C 1 - allowing for a controlled perturbative expansion at 
low momenta, k < A (see Fig. 1). Conceptually, our approach therefore offers a well-defined 
and controlled treatment of the effects of short-distance non-linearities on the long-wavelength 
universe. This is to be contrasted with the failure of many cosmological perturbation theory 
techniques to include the effects of gravitational non-linearities [20]. 

conventional perturbation theory ^ 



effective theory ^ 



NL 



Figure 1: Hierarchy of scales and perturbative expansions. In the effective theory loop integrals only 
contain modes with k < A, while conventional perturbation theory contains modes with k ~ k^ where 
the perturbative expansion is known to break down (see §6j. 



4 We point out that the effective viscosity of the fluid leads to a damping in the non-linear evolution of density 
fluctuations 8i and suggest that this intuitively explains the non-linear broadening of the peak of baryon acoustic 
oscillations (see §6). 
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Outline: The outline of the paper is as follows: In Section 2 give a bottom-up construction 
of the effective stress-energy tensor using the Newtonian approximation on subhorizon scales. 
Alternatively, in Section 3 we explain our basic perturbative approach within General Relativ- 
ity and present second-order results for metric and matter perturbations in Poisson gauge. In 
Section 4 we derive the gravitational stress-energy pseudo-tensor and define an effective fluid 
by taking its long- wavelength limit. We show that on very large scales the fluid behaves as 
an isotropic fluid and gravitational non-linearities only renormalize the background density and 
pressure. Furthermore, we prove that virial scales decouple completely from the long-wavelength 
theory. In Section 5 we show that on scales comparable to and smaller than the horizon scale 
the fluid anisotropic stress is non-negligible and important for the evolution of perturbations. 
We characterize the properties of this imperfect fluid in detail. In Section 6 we suggest possible 
applications of our effective theory. Finally, we present our conclusions in Section 7. 

A number of appendices contain technical details: Appendix A gives an alternative derivation 
of the effective fluid properties starting from the Newtonian conservation equations. This deriva- 
tion provides considerable intuition for the physical origin of the effective fluid. Appendix B 
presents more details of second-order cosmological perturbation theory. We collect results in 
Poisson gauge, cite second-order gauge transformations and discuss the long-wavelength limit of 
the spacetime. In Appendix C we review key elements of dissipative fluid dynamics that are 
used in Section 5. Finally, in Appendix D we illustrate some of our ideas by presenting example 
calculations in perturbation theory. 

Short Proofs and Examples are separated from the main text by horizontal lines. These parts can 
be omitted without loss of continuity, but they often illuminate the underlying physics. 

Our work builds on a long history of studies in non-linear perturbation theory [23, 24, 25, 26, 
27, 28, 29, 30, 31, 2], numerical simulations [32, 33, 34, 35], and effective field theory [5, 6, 8, 37]. 
Related ideas have appeared recently in Refs. [38, 39, 40, 41, 43, 44, 45]. 

Notation and conventions: Except for the Newtonian analysis of §2, we will work exclusively 
with conformal time r\ and in units where the speed of light is set to unity, c = 1. Greek 
indices [i, v = 0,1,2,3 are used for four-dimensional spacetime coordinates, while Latin indices 
i,j = 1,2,3 are reserved for spatial coordinates. We will use overdots to indicate derivatives 
with respect to conformal time. We use commas to denote partial derivaties and semi-colons for 
covariant derivatives, i.e. 



■) = »/.(•■■) 



and 



(...) 




(1.7) 



Our Fourier convention will be 




(1.8) 



where we used the notation 




(1.9) 
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Reality of <p(x) demands that <$>* k = The product of two functions in real space is a convo- 
lution in Fourier space 

f e- ik - x <t>{x)^{x)= [ <f> q i> k - q . (1.10) 

J x J q 

We define two forms of the power spectrum, 

{Mv)Mv)) = (2TT) 3 5 D (k + k')P4k, v ), (1.11) 

and 

^l(k,v) = ^P<t>(k,v), (1-12) 

such that (</> 2 ) = J din k A 2 (fc). The time-evolution of the gravitational potential on subhorizon 
scales during the radiation era leads to a momentum-dependent transfer function T^k, rj), where 

Mv) = T ( p(k,r l ) ( t> k (0). (1.13) 

We take the initial conditions to be scale-invariant, with an amplitude fixed by cosmic microwave 
background observations [46], A|(&,0) « 10 -9 . 



2 Bottom-Up Construction of the Effective Stress-Energy 

As we stressed in the previous section, the large hierarchy between the non-linear scales and the 
horizon size, and the smallness of typical velocities and gravitational potentials for the former, 
make the problem of 'backreaction' amenable to a Newtonian analysis. To make use of the fact 
that the expansion of the universe is negligible at subhorizon scales, it is particularly convenient 
to use coordinates for which locally around any given comoving point the unperturbed FRW 
geometry is manifestly a small perturbation of Minkowski space [47, 48] . That this is possible is 
a straightforward consequence of the equivalence principle. However, what is not that straight- 
forward is that we can do so for an entire comoving geodesic — i.e. the metric is approximately 
flat at subhorizon distances from the geodesic, for all times. 5 This is important because we do 
not want to limit the validity of our approximation to times shorter than H^ 1 (since H^ 1 is the 
time-scale we will eventually be interested in). Taking x = as the preferred comoving point, 
the background FRW metric locally becomes [48] 



d4 RW - 



l-(H + H 2 )x 2 ] dt 2 + [l - \h 2 x 2 ] dx 2 = (^ + h™) dx^dx" , (2.1) 



where x and t are suitably defined 'physical' (as opposed to comoving) coordinates, and H 
and H are evaluated at t. The actual relation between (t, x) and the standard comoving FRW 
coordinates (t c , x c ) is 

t c = t-\H(t)x 2 , x c =^[l + \H 2 (t)x 2 ]. (2.2) 



5 Local coordinates with these properties are called Fermi coordinates [47]. 
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In these coordinates the Hubble flow corresponds to v ~ Hx. Corrections to the above metric 
are suppressed by further powers of ~ H 2 x 2 . We thus see that when studying a non-relativistic 
system of subhorizon size, we can encode the Hubble expansion into a Newtonian potential 6 

$frw = -\{H + H 2 )x 2 (2.3) 

and a background velocity field 

■ufrw = Hx . (2.4) 

Perturbations of the background homogeneity appear as further contributions to the total $ and 
v fields. 

For simplicity, let us imagine the universe to be filled with a pressureless fluid. 7 In the New- 
tonian approximation the dynamics of the fluid coupled to gravity is governed by the continuity 
and Euler equations, supplemented with Poisson's equation for the Newtonian potential, 

Pm + V- (p m v) = 0, (2.5) 
v + (v-V)v = -V$, (2.6) 
V 2 $ = 4vrG / 9 m . (2.7) 

For instance, the background fields (2.3) and (2.4) obey these equations exactly, for homogeneous 
p m , provided that H 2 and H satisfy the usual Friedmann equations. We stress that here p m is the 
fluid's Newtonian mass density. We use this notation to be consistent with the rest of the paper, 
where we stick to the standard relativistic convention of reserving p for the fluid's energy density 
as measured in its local rest-frame. Our aim is to construct a symmetric stress-energy tensor 
t^ u that is conserved by virtue of the above equations of motion. We are guided by knowing the 
conserved total Newtonian energy and momentum of the system 

J d 3 x [p m + X -p m v 2 + \pm®\ = j d 3 x r 00 , (2.8) 

P i = J (Pxprrrf = J d 3 XT 0i . (2.9) 

Of course, these integral expressions make sense only for finite systems — for infinite systems like 
the universe, they are divergent. We use them here just as a guide to formulate an ansatz for the 
local stress-energy tensor. Once we explicitly check that the stress-energy tensor thus guessed is 
conserved, we have no need to define global quantities like the total energy and momentum. Let 
us thus start by postulating 

T° l = pmv' (2.10) 

and by imposing the i-component of the conservation equation 

= fl/ = Pm v l + p m v % + (2.11) 

= -dj ( Pm »V) - p m di$ + dj^ , (2.12) 

6 For subhorizon relativistic systems we would also need the second Newtonian potential $frw = \H 2 x 2 . 
In this section we are using the standard convention for the metric in the Newtonian approximation: ds 2 = 
-(1 + 2$)d£ 2 + (1 - 2*)dx 2 . 

7 We discuss the limitations of this approximation in Appendix A. 



E = 



10 



where we used the continuity and Euler equations. If we are able the rewrite the second term 
also as a divergence, we have an expression for r y . This is straightforward to do upon using 
Poisson's equation 



1 



AttG 

We thus have a natural candidate for r y : 



;V 2 <^ = -^dj^dj* - ^-(V<&) 2 



A3 = - 



p m v l v 3 + 



1 



8ttG L 



2d i $d <7 -$ 



(2.13) 



(2.14) 



For r 00 instead, we have several natural choices. This is because, as usual, the gravitational 
contribution to the energy density in (2.8) can be written as p m & or as (V<£) 2 or as any linear 
combination thereof — these possibilities all coincide upon using Poisson's equation and integra- 
tion by parts. However, this ambiguity corresponds precisely to higher-derivative terms of the 
form (1.1). Indeed 

(V$) 2 -> -4irGp m $ + \v 2 § 2 . (2.15) 



The associated £ tensor — cf. Eqn. (1.1) — is thus 



(2.16) 



and zero for all other components not related to these by symmetry. We find it more convenient to 
use the (V<E>) 2 representation for the gravitational potential energy. Our ansatz for r 00 therefore 
is 



-00 



Pm + ^PmV 



8vrG 



We now impose that this obeys the 0-component of stress-energy conservation 

= 8^° = d T 00 + d iT i0 . 



(2.17) 



(2.18) 



We do not assume that r*° is the same as r° l defined in (2.10). As we will see in a moment this 
is an interesting point. Taking the time-derivative of (2.17) and using repeatedly the continuity, 
Euler, and Poisson equations, we get 



dor 00 = -di 



AttG 



This is consistent with the local conservation law for 

1 



jo 



= PmV % ( 



1 + + $ + 



(2.19) 



(2.20) 



Now, as we have just shown, the stress-energy tensor defined by Eqns. (2.10), (2.14), (2.17), 
and (2.20) is exactly conserved for our Netwonian system and yields the usual total energy 
and momentum of Netwonian mechanics. However r°* 7^ r l °, so before using t^ u as a source 
for long-wavelength gravitational fields, we should make it symmetric. Symmetrization of the 
stress-tensor is always possible (see e.g. Ref. [17]), but it requires Lorentz-invariance because it 
explicitly involves the use of the Lorentz algebra generators. This means that we cannot make 
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the stress-energy tensor symmetric within our non-relativistic approximation. This is hardly a 
problem however, since the mismatch between r*° and r° l is negligible at non-relativistic speeds 
and gravitational fields. This is evident for the terms in parentheses in Eqn. (2.20), whereas for 
the last term we have 

4^® 9i ® ~ ~4^G^ Vdidj ® ~ PmV^^pmV. (2.21) 

This means that in our non-relativistic Newtonian approximation we can use our t^ v with r l ° set 
equal to r° l = p m v l , as the post-Newtonian source for gravity. Notice in passing that we never 
required inhomogeneities in p m to be small. In fact, they can be huge, as long as velocities and 
gravitational fields are non-relativistic. 

A few remarks are in order: For a Newtonian system governed by Eqns. (2.5)— (2.7) there is 
another locally conserved and symmetric tensor f^", 

f 00 = p m , f 0i = r i0 = p m v l , (2.22) 

and = T l i as given in Eqn. (2.14). The 0-component of its conservation equation is precisely 
the continuity equation, which expresses the conservation of mass. Indeed the global charge 
associated with f 00 is the total mass of the system 

73 



M = j <Txpm . (2.23) 

Now, in Newtonian physics mass and energy (kinetic plus potential) are separately conserved. 
However, while the latter is conserved thanks to the dynamics, the conservation of the former is 
a trivial kinematical constraint — in the sense that the mass of a system is part of its definition 
rather than being a dynamical quantity whose value depends on the system's configuration. 8 We 
know that upon inclusion of relativistic corrections this degeneracy is gone, and only one form of 
energy is conserved, namely that which in the Newtonian limit reduces to (2.8). This guarantees 
that the stress-energy tensor that sources the gravitational field is indeed t^ v rather than f ^ (or 
any linear combination of the two). This conclusion will be manifest in the general-relativistic 
treatment of §3. 

Next, we want to consider the contribution of the FRW background fields (2.3) and (2.4) to 
the effective stress-energy tensor. First, in the absence of fluctuations we schematically have 

Tfrw ~ Pm (1 + H 2 x 2 ) , 4i w ~ p m Hx , 4 J RW ~ Pm H 2 x 2 . (2.24) 

The explicit x-dependence is clearly related to our choice of coordinates, which is not particularly 
convenient to address questions beyond leading order in Hx. Of course, the FRW geometry is 
invariant under properly defined translations, and so should be any suitably defined physical 
quantity, but this is not manifest at all in our coordinate system. Likewise, when we have 
fluctuations in $ and v, 

$ = <&frw + , v = i>frw + Sv , (2.25) 



8 The formal quantum-mechanical counterpart of this statement is that the mass appears as a central charge in 
the algebra of Galilean transformations and translations [17]. As a consequence, states of different mass belong to 
different super-selection sectors of the theory. 
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the mixed background- fluctuation terms in are weighed by Hx. This yields a stress-energy 
tensor for the fluctuations that depends explicitly on position. Given the background physical 
homogeneity, we would like to get rid of this fake explicit position-dependence and be left just with 
the physical one implicitly contained in the fluctuation fields. We want to stress that as long as 
we stay sufficiently close to the origin so that special- and general-relativistic effects are negligible, 
our expression for t^ u is certainly the correct one in the coordinates we have been using, even 
in the presence of an FRW background. Only, such coordinates hide the nice symmetries of the 
cosmological background and are thus not suitable for cosmological applications. Of course, in 
the rest of the paper we will use standard FRW coordinates — so how do we translate our t^ u to 
FRW comoving coordinates? One possibility is to perform the actual coordinate transformation, 
but the gravitational part of the stress-energy tensor, infamously, does not transform as a tensor 
under generic coordinate transformations. Indeed, it can always be made to vanish at any given 
point by choosing a locally inertial frame [19]. Another possibility, which we will adopt, is to 
evaluate our t^ u at x = 0. That is, given the homogeneity of the background, any comoving 
point P is as good as any other to play the role of x = 0. Therefore, if we are interested in the 
components of t^ v at P in FRW coordinates we can proceed as follows: We identify the origin 
of our Netwonian coordinates (2.1) with P. We compute t^ u at x = according to Eqns. (2.10), 
(2.14), and (2.17). The background fields (2.3) and (2.4), given their ^-dependence, give no 
contribution at x = 0. That is, at the origin only the fluctuations in <3? and v contribute to 
t^ u . Finally, we perform the coordinate change (2.2) which brings us back to comoving FRW 
coordinates, and which at the origin is exceedingly simple. There is a factor of a(t) for each spatial 
index, and more importantly, the problematic non-tensor part of the transformation law for t^ u 
is gone, because the point x = is at rest for all times in both coordinate systems — its state 
of motion is untouched by the coordinate change. In conclusion, when we go back to standard 
FRW coordinates, we can use our expressions for t^ v directly in terms of the Newtonian-gauge 
potential and peculiar velocity, modulo obvious factors of a(t). 

Finally, it is instructive to present the derivation of our effective stress-tensor in the 
Newtonian context in yet another way. As we mentioned in the Introduction, we will later define 
the effective theory for long-wavelength fluctuations by smoothing the stress-energy tensor t^ u 
on a scale A and declaring that long- wavelength gravitational fields are coupled to it. It is 
particularly illuminating to see how t^ v arises in we perform the smoothing immediately at the 
level of the Euler and Poisson equations (2.5) and (2.6). We apply a filter on scales of order A -1 
to the Euler equation 



and split the fields into short-wavelength and long- wavelength fluctuations X = X^+X s . Straight- 
forward algebra then shows (see Appendix A) that the Euler equation can be recast in the fol- 
lowing way 




(2.26) 



We define smoothed quantities of all fields X = {p m , <&, p m v} as 




(2.27) 



(2.28) 
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where 

[nj] 3 = [ Pm v!v°] A + — [2d^ s d^ s - ^(V$ s ) 2 ] A . (2.29) 

We see that the long- wavelength fluctuations obey an Euler equation in which the stress tensor Tij 
receives contributions from the short-wavelength fluctuations. Eqn. (2.29) is exactly of the form 
of the effective derived in Eqn. (2.14), and it shows explicitly how the effective long- wavelength 
fluid is different from the pressureless fluid we started with in the continuity and Euler equations 
(2.5) and (2.6). 

The skeptical reader will find comfort in the analysis of the next section, where we derive 
the effective stress-energy tensor directly in FRW by more traditional means. 9 Those already 
convinced that our t^ v is the correct post-Newtonian source for gravity, can instead skip directly 
to §4 without substantial loss of continuity. For these readers, the starting point should be 
Eqns. (3.25) and (3.26). There, t^ v is expressed in terms of p, the energy density as measured 
by an observer comoving with the fluid. This differs from our p m by terms of order pv 2 and /9$. 
Also, the t 00 given in Eqn. (3.25) differs from that of this section by a different choice for the 
p<& vs. (V<I>) 2 ambiguity in the potential energy. As we argued different choices correspond to 
different higher-derivative terms of the form (1.1), which are irrelevant at long wavelengths. 

3 Perturbation Theory and Velocity Expansion 

We now turn to the more standard approach to cosmological perturbation theory, where one 
expands in perturbations around the FRW metric. However, contrary to the standard way of 
organizing the expansion in fluctuations, we will show that perturbation theory may be formu- 
lated systematically as an expansion in perturbed velocities, in analogy with the post-Newtonian 
approach to General Relativity. This has certain advantages when trying to capture non-linear 
perturbations. In §3.1 we describe the general philosophy of this approach. We introduce the 
concept of the effective stress-energy tensor arising from short- wavelength perturbations. The 
treatment will be schematic with details postponed to the following sections. In §3.2 we give the 
equations of second-order perturbation theory in Poisson gauge (see also Appendix B). This is 
used in §3.3 to derive the effective stress-energy tensor. 

3.1 Generalities 

Consider the perturbed spacetime 10 

ds 2 = a 2 {r]) [_ e 2* d?? 2 + g -2$ d;r 2] ; {31) 

In comparing the results of §2 and §3 the reader should be warned that in §3 we will switch the convention for 
the metric fluctuations $ and ^> (see Eqn. (3.1)) in order to follow the standard convention for Newtonian gauge 
(to be distinguished from the convention used in the Newtonian approximation in §2). The precise correspondence 
is: <5$§2 = *&§3 , <5^§2 = $§3. In comparing the results for the stress tensor r M „ at second order in $ and ^ one 
should also use that in our approximations the first-order perturbations are equal: = ^§3'- 

10 For the purposes of this introductory section we neglect vector and tensor perturbations and we don't worry 
about technical details like gauge-fixing. These deficiencies will be cured momentarily. 
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sourced by matter perturbations Sp, Sp and V{. In standard cosmological perturbation theory the 
Einstein equations are expanded in metric perturbations and matter perturbations. Perturbation 
theory is then said to break down when the magnitude of density fluctuations Sp(r], x) becomes 
comparable to the background value pirf). However, this gives up too early. It is well-known 
that metric perturbations and matter velocities remain small even if the density perturbations 
become large, dp > p. Indeed, the density contrast S = Sp/p is ~ 10 2 , 10 10 and 10 30 on cluster 
scales, galaxy scales and solar system scales, respectively, while the spacetime on these scales is 
still well-approximated by a homogeneous spacetime with small perturbations $ < 1CP 5 . 

The following alternative to conventional cosmological perturbation theory therefore suggests 
itself: the Einstein equations are expanded in metric perturbations <£, ^ and matter velocities 
v % , but not in density perturbations Sp. In perturbation theory the linear peculiar velocities are 
related to the Newtonian potential and the density contrast as follows 

v 2 ~ $<5. (3.2) 

When both <3? and S are much smaller than unity, v 2 is therefore considered a second-order 
quantity. However, when entering the quasi- nonlinear regime, 5nl ~ 1, the fluid kinetic energy 
(per unit mass) effectively becomes a 'first-order' contribution, v 2 ~ <£. Similarly, for virialized 
objects, (5 v ; r ~ 10 2 , one also finds that kinetic and gravitational energies are of the same order, 
v 2 ~ <E>. 

Furthermore, on small scales, gradients of the gravitational potential can change the power- 
counting of standard perturbation theory. For instance, quadratic terms with two spatial deriva- 
tives satisfy the following scaling relations 

(V<I>) 2 $V 2 $ 

~ ~ ~ v 2 , (3.3) 

where H = a/ a is the comoving Hubble parameter and we have used the Poisson equation (for a 
pure matter universe), V 2< 5 = |% 2 5. For small-scale perturbations the terms (V^) 2 and <I>V 2( I> 
are therefore enhanced relative to terms without at least two gradients 11 , like % 2< 3? 2 . For a (scale- 
invariant) spectrum of perturbations with different wavelengths, terms with the maximal number 
of gradients will hence give the leading effect. We will capture this by taking v 2 as our expansion 
parameter. At the non- linear scale <I> 2 is order v A while each gradient (rendered dimensionless 
by the Hubble scale %) reduces the order in v by one. This is reminisent of the post-Newtonian 
approach to General Relativity. 

Of course, our statement about expanding in velocities is not gauge-invariant: one could even 
choose a comoving gauge where velocities are zero, in which case clearly our expansion would not 
be well-defined. However, this is not a source of worry for the following reason: the fluctuations 
that become non-linear are those related to the matter degrees of freedom, while the metric 
fluctuations are not large even on very short scales (except close to black holes 12 ). Of course, one 



11 In General Relativity, the equations of motion have at most two derivatives. This implies that the perturbation 
expansion naturally truncates at second order and third-order terms are suppressed by an extra power of $, 
e.g. $ 2 V 2 $ < $V 2 $. 

12 The case of order one metric fluctuations near black holes can also be taken into account quite straightforwardly. 
We have briefly commented on this in the Introduction. 
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can choose a gauge in which the matter fluctuations vanish, and all the fluctuations are 'eaten' by 
the metric. In this gauge metric fluctuations will be large. In this case, the large fluctuations of 
the metric simply are related to a contorted slicing of the spacetime, which is, in fact, never very 
different from FRW. This particular fact implies that there exist gauges where metric fluctuations 
are small and matter fluctuations can be large. This is the case for Newtonian gauge, or for any 
gauge where the matter fluctuations are not made explicitly small by the gauge condition. Our 
expansion scheme applies in those gauges. A gauge where instead metric fluctuations are as 
large as the matter fluctuations is the so called £-gauge, or unitary gauge. In this case, matter 
fluctuations are actually taken to be zero, and our expansion is not applicable. 

The reason why there are universes that at subhorizon distances feature large matter fluctu- 
ations but small metric perturbations can be stated in a language that is familiar from particle 
physics where is goes by the name of the equivalence theorem for massive vector fields. A non- 
trivial matter background density spontaneously breaks some of the spacetime symmetries, which 
are gauged in General Relativity. As a consequence the gravitational field gets a mass of order 
Hubble. Matter fluctuations can be considered as Goldstone bosons of these symmetries, which 
according to the equivalence theorem, are the most strongly-coupled degrees of freedom at scales 
smaller than the graviton inverse Compton wavelength, that is at subhorizon scales. This is 
because their self-couplings involve more derivatives than those they have with the metric fluc- 
tuations or than the metric fluctuations' self-couplings. As a consequence, on small distances 
matter fluctuations are more strongly coupled, and therefore non-linear, than the metric ones. 
This hierarchy is completely obscure in the gauge in which matter fluctuations are set to zero — 
in which they are 'eaten' by the metric — which is called unitary gauge in the particle physics 
context. This logic is expanded upon in Ref. [36, 37]. 

As we just remarked, our approach of expanding in velocities is valid for all gauges in which 
there is an hierarchy between matter fluctuations and metric fluctuations. One can consistently 
change gauges within this set (see Appendix B). Our approach will lead us to neglect terms 
quadratic in the metric perturbations of order $ 2 ~ 1CP 10 . Because of this, our approach is not 
gauge-invariant at this level. Furthermore, at this order it even becomes unclear how to split the 
metric into long- wavelength and short-wavelength fluctuations (see §6.4), as we did in Eqn. (2.27) 
and will do in the rest of the paper. Luckily, there are many non-linear effects in the matter 
sector whose size is larger than 1CT 10 . It is those effects that can be addressed unambiguously. 
It is unclear to us if effects of order <E> 2 ~ 10 -10 will ever be relevant from the observational 
point of view in the context of LSS surveys. In any case, describing those effects would require 
a qualitative improvement of our treatment that goes beyond the scope of the current paper. 

We are now ready to return to the main objective in this work, which is to develop an effective 
description for the evolution of long- wavelength perturbations in the presence of short- wavelength 
non-linearities valid for effects larger than order $ 2 ~ 10~ 10 . To achieve this, we decompose the 
Einstein tensor into a homogeneous background (denoted by overbars) and terms that are linear 
(L) and non-linear (NL) in the perturbations 5X(rj, x) = X(rj, x) — X(r]). The Einstein equations 
can then be written as 

[X] + (G M ,) L [SX] + (G M ,) NL [5X 2 ] = 8itG . (3.4) 
The background equations, G^ u = 8irGT^ u , and the linearized Einstein equations, (G F ) L = 
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8irG (T M j,) L , are defined in the standard way. The non-linear Einstein equations can be written 
in a form that is very similar to the linear equations, 

(G^) l = 8wG(t^-T^), (3.5) 



where we defined the effective stress-energy pseudo-tensor 

r, v = T, u - • (3-6) 



In the following we compute the effective stress-energy pseudo-tensor associated with the short- 
wavelength non-linearities at quadratic order and as an expansion in the matter velocities. 



3.2 Second-Order Perturbation Theory 

We now start to define things more precisely. We begin by collecting standard results in second- 
order perturbation theory in Poisson gauge, a generalization of Newtonian gauge to second order. 

Metric perturbations. The metric in Poisson gauge is given by 

ds 2 = a 2 {r,) [-e^dr/ 2 + 2^dxMr/ + (e~ 2 %- + X ijW l dx j ] , (3.7) 

with xa = — i.e. the trace of gij is absorbed into $ — and 

= Xij,i = • (3.8) 

The gauge condition (3.8) eliminates one scalar degree of freedom from go* an d one scalar and 
one transverse vector degree of freedom from gij. Thus, Ui is a transverse vector, while Xij is a 
transverse-traceless tensor. The metric determinant is y/^g = e* -3 *^ up to second order in per- 
turbations. Expressions for the Christoffel symbols, Riemann and Einstein tensors corresponding 
to the spacetime (3.7) are given in Appendix B (see also Refs. [50, 51, 52]). 
Perturbations are formally split into first-order and second-order terms 

* = * (1) + ^ (2) , * = fcW + ±*( 2 > , c* = W < 2 >, x^ = xg } . (3.9) 

Here, we have ignored first-order vector and tensor perturbations, i.e. = x[f = 0. Many 
models of inflation indeed produce initial conditions where vector modes are zero and tensor 
modes are negligibly small [53, 54]. Furthermore, at first order, subhorizon vectors and tensors 
decay as the universe expands. The absence of first-order vectors and tensors has the following 
important consequence: second-order scalar modes are only sourced by first-order scalars. In 
fact, although scalars, vectors and tensors mix beyond linear order in perturbation theory, the 
second-order parts of scalars, vectors and tensors do not mix, e.g. second-order scalars are only 
sourced by quadratic contributions from the first-order vectors and tensors. 
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The non-linear parts of the components of the Einstein tensor are (see Appendix B), 



\NL 



V(G U 
a 2 ■ 



,NL 



12H 2 * 2 + 12H$y + 3<£ 2 - <S> ik <S>,k + 4$$ fefc , (3.10) 

2$[$ + H*], i -6* i , (3.11) 

-4(H 2 + 2ft)* 2 - 2$tf - 3<I» 2 - 4-H(2$ + ^)^ - 4*$ 

+ * fc * fc -2*[*-*] )fcfc ](5j 
+ 2<l> <l< - * . u + <l>,<l>., - * ti *j - *;<!>.; - . (3.12) 



To define matter perturbations we introduce the timelike velocity four- 



Matter perturbations. 
vector 

,-/,'/' 

(3.13) 

where r is the proper time comoving with the fluid, so that g^ v u^u u = —1. We define the 
tensor 7^ = g^ u + u^u u which projects tensors orthogonal to the four-velocity into the fluid's 
instantaneous rest space at each event. The energy-momentum tensor of a general (imperfect) 
fluid decomposes, at leading order in derivatives, into irreducible parts as 



Tfus = pu^u v + (p- C#h^ + > 



(3.14) 



i lu u fJ 'u u is the matter energy density, p is the isotropic pressure, ( is the bulk 



where p = T f 

viscosity, 9 = diUi is the velocity divergence and T,^ = J^jJ}T a p is the symmetric and trace- 
free anisotropic stress tensor. 13 and ( vanish for a perfect fluid, ( = = 0. It holds that 
p - (8 = \T^ U 

For cold dark matter we furthermore impose that the background pressure vanishes. Our 
starting point is therefore the stress-energy tensor, 15 



T% = pv! 1 u v . 



(3.15) 



To second order in metric and three- velocity perturbations the components of the fluid four- 
velocity are 

u° = a- 1 e- 9 'y(v), u^a^eV, (3.16) 



where we defined the Lorentz factor 



7(1;) = 



1+ 2 V 



(3.17) 



13 Here we use the notation = 7(^7^*a^ ~ §7 Q ^*a/37*»' and = \ {t^ v + t vll ). 

14 It is worth pointing out that the viscosity is a higher-derivative term that locally in space and time is indis- 
tinguishable from the pressure. This explains why we cannot isolate the pressure through a tensorial contraction. 

15 On very small scales the fluid description is expected to break down. In that limit we will treat the dark matter 
as a collection of point particles that interact only via gravity [18] (see Appendix B). For notational simplicity we 
will often use the perfect fluid description of the dark matter, i. e. the continuum limit of the particle model and the 
S/aK — P = limit of Eqn. (3.14). However, all our results are easily generalized to the background of particles and 
sometimes we will find it more convenient to think of the dark matter in terms of particles (e.g. when discussing 
the effects of tidal forces from long-wavelength perturbations in §5.2). When relating our approach to cosmological 
iV-body simulations of dark matter (see §6) the effective stress-energy is that of non-relativistic point particles. 
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At order v 2 the components of the stress-energy tensor during matter-domination therefore are 
T° = - 1 2 p = -p{l + v 2 ), r i = -e*+V, T)= P v\. (3.18) 
We emphasize again that we do not expand in density fluctuations. We may still write 

3W 2 

p(r),x) = p(rj)[l + 6{T},x)], where p = gjrGa2 , (3-19) 

but we do not assume that 5 < 1 (unless we are specifically referring to long-wavelength density 
perturbations). 

Einstein equations. The metric and matter perturbations are coupled to each other by the 
Einstein equations. As we have argued above, the non-linear Einstein equations can be written 
in a form that is very similar to the linear equations if we replace the stress-energy energy tensor 
Tt by the stress-energy pseudo-tensor t%: 

(G^) L = SttGK - Tt) , (3.20) 

or explicitly for scalar fluctuations, 

V 2 <5>-3n{$ + ny) = -4vrGa 2 (r° - f° ) , (3.21) 
+ = 4vrGaV , (3.22) 

§ + U{2$> + i>) + {U 2 + 2U)^ - \v 2 {<S>-^>) = ^^-(t\-t\), (3.23) 



didj 



($-%-^iV 2 ($-$) 



87rGa 2 didj 



(3.24) 



This defines the dynamics of long-wavelength scalar fluctuations sourced by products of short- 
wavelength fluctuations. All gravitational non-linearities have been moved to the r.h.s. of the 
Einstein equations to define the effective stress-energy tensor t%. The tensor rt is conserved by 
virtue of the linearized Bianchi identity — see Eqn. (B.69)). 

3.3 Effective Stress-Energy 

From the expressions for the second-order Einstein tensor, Eqns. (3.10)— (3.12), we read off the 
components of the stress-energy pseudo-tensor 

T o = -p{l + v k v k ) 87T Ga 2 ' ^ ' 

r, = „ ^ . (3.26) 

Here we have used the fact that = = 0. Corrections to r°j are of order w 3 and hence 
not shown. We define the trace r = r % i and the traceless part f l - = — |5jr^.. The divergence 
of Eqn. (3.26) is 

ViT-V = V;(p^) + ~ Vi(puiUj) + . (3.27) 

Eqns. (3.25) and (3.26) form the basis of our exploration of gravitational non-linearities as an 
effective fluid. 
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4 Non-Linear Gravity as an Effective Fluid 



Having derived the effective stress-tensor associated with non-linear fluctuations in two different 
ways (cf. Sections 2, 3 and Appendix A), we are now in the position to discuss in detail the effects 
of short-wavelength fluctuations on the long-wavelength universe. 

In §4.1 we explain the procedure of integrating out short-wavelength fluctuations to arrive at 
the long-wavelength limit of the effective stress-energy pseudo-tensor t^. In §4.2 we digress to 
explain why higher-order moments in the Boltzmann hierarchy are suppressed for dark matter 
particles in a FRW universe. This truncation of the Boltzmann hierarchy allows us to describe the 
system at long wavelengths as a fluid with density, pressure and anisotropic stress. We describe 
the properties of t^ u in terms of an effective fluid in §4.3. In §4.4 we show that anisotropic stress 
is negligible on very large scales and that integrating out small scales is therefore captured simply 
by a renormalization of the density and pressure of the background. However, on scales smaller 
than the horizon, anisotropic stress becomes important for the evolution of long-wavelength 
perturbations. We will describe the parameters of this imperfect fluid in the next section, §5. 

4.1 Integrating out Small Scales 

The leading contributions to the non-linear gravitational source terms, Eqns. (3.25) and (3.26), 
contain two spatial derivatives. For scale-invariant perturbations in the Newtonian potential 
the effects are hence dominated by small scales with characteristic momentum q*. Since we are 
interested in the theory at scales k much larger than the scale of non-linearities q*, we define 
an effective long-wavelength theory by 'integrating out' short-wavelength modes below a scale 
A <C q+. Here, integrating out short-wavelength fluctuations amounts to smoothing the equations 
of motion and taking expectation values of the short- wavelength modes in the presence of long- 
wavelength perturbations, so that one is left with equations in terms of only the long-wavelength 
modes. 16 In real space, the smoothing of perturbations corresponds to a convolution of all fields 
X = {p, <!>, ^, pv} with a window function W A , 

X e = [X} A (x) = J d 3 x' W A (\x - x'\)X(x') . (4.1) 

Essentially, Eqn. (4.1) amounts to averaging the fields over domains of size A -1 . Note that we 
smooth the momentum density, j = pv, rather than the velocity v. From a dynamical point of 
view, it is more natural to average the momentum rather than the velocity, since the former is 
an additive quantity for a system of particles. We then split all fields into long and short modes, 

X = X e + X s . (4.2) 

Next, we consider the smoothing of the Einstein equations (3.21) - (3.24), 

J d 3 x' W A (\x -x'\)- Gl u (x') = J d 3 x' W A (\x -x'\)-{ - } (4.3) 

= Wa-T f . (4.4) 

16 Taking the expectation value of the short modes isn't strictly necessary. The backreaction of the short- 
wavelength modes on the long-wavelength modes could also be dealt with also on a realization by realization 
basis. 
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Being linear in metric fluctuations and the l.h.s. of Eqn. (4.3) extracts the long-wavelength 
metric perturbations and To compute the non-linear source terms for 3>£ and \E^, we 
substitute the fields in Eqn. (4.2) into the r.h.s. of Eqn. (4.3), using Eqns. (3.25) and (3.26) for 
the effective stress-energy tensor t^ u . We find, 

[v]A = [vfl[v] s + [r/, (4.5) 

where [t^Y depends only on long-wavelength modes, Xe, and [t^ u ] s is quadratic in short- 
wavelength modes, X s . Mixed terms with one long- wavelength mode and one short-wavelength 
mode vanish up to higher-derivative terms which are suppressed by powers of k/A. Collectively 
we denote all higher-derivative corrections by [t^]® 2 . Typically, these higher-derivative terms 
can be dropped unless the aim is very high precision. Non-linear metric contributions in [t^Y 
may be moved back to the l.h.s. of the Einstein equations to define the non-linear completion 
of the Einstein tensor for long-wavelength modes, [G^Y = [G^ V Y + [^JSt- These modes are 
sourced by [t^] 3 . 

After a computation completely analogous to that performed explicitly in Appendix A (see 
also the Example below) we find: 



• long-wavelength modes 



r om n , k tk^k-^^kk 

W = -Mi + 44) - 87TGa 2 > ( 4 - 6 ) 



JJ ' i 3 ^Ga2 

• short-wavelength non-linearities 



Kr—M-^^^, (4.8) 

^. lP ^m^mu. (4 . 9) 



• higher- derivative terms 



[T ° ] = ~ Pe ^^ SvGa? ■ A^ ' (410) 

[T i ] - pe ^^ + SvrGa* • A* " (4U) 



Example: 

The derivation of Eqns. (4.5) - (4.11) is tedious, but instructive. In order to provide some intuition for 
the physical origin of the different terms we now give a sample computation, which applies to all bilinear 
terms in t m „. All other terms (i.e. the trilinear ones <~ pvv) work in a similar way and are shown in 
Appendix A. 
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Consider smoothing of a bilinear quantity fg, 

[fg]A= [ W A f(x')g(x'). (4.12) 
We split the fields /, g into long- wavelength modes ft, ge and short- wavelength modes f s , g s . We get 

[fgh = [fmh + [f s g s h + [fegsh + [f s geh ■ (4.13) 

The second term is already in the short-short form of Eqns. (4.8) and (4.9). We now massage the other 
terms. In the first term, given their mild dependence on x' , we can expand the long modes fi(x') and 
ge(x') in a Taylor series about x: 

lftgeh(x) = f W A f e {x')g e (ae?) (4.14) 

J x' 

= hat + ^2 i^h ■ V» + \feV 2 g t + \g,V 2 fi) +... , (4.15) 

where the dots stand for even higher derivative terms. We used the normalization condition for the window 
function, J W A — 1, and the characterization of 1/A 2 as the real-space variance of W A : 



L 



W A (x-x'y(x-x'y = ^S ij . (4.16) 



In fact, we take this as the definition of what we mean by the smoothing scale A 1 . This way the expansion 
(4.14) docs not depend on the specific window function we adopt (we are of course assuming isotropy of 
Wa). 

To simplify the mixed long-short terms in Eqn. (4.13) we have to work a little harder. For instance, 
for the third term we rewrite 

[fegsh = [figh - [ftgth ■ ( 4 - 17 ) 

Here, the second term we already dealt with, so we have to simplify the first term. However, to proceed 
we need further assumptions on W\ - i.e. the result will depend explicitly on the form of the window 
function. A particularly convenient choice is a Gaussian, because then we have useful identities: 

d f W A = -djW A = A 2 {x - x') j W K , (4.18) 
di>d f W A = didjW A = -A 2 SijW A + K A {x - x') l (x - x') j W A , . (4.19) 

We can then expand fg in a Taylor series about x like before and get 

[fegh = fm-Vfi-[{x-x')g(x l )h + \v i V ,f i -[{x-x l )\x-x l y g(x')] A + .-. (4.20) 

= fegi + (V// • Vg e + \g f y 2 f t ) + . . . (4.21) 
In conclusion, our smoothed bilinear is 

[fgh - ftgt + [fegsh + ^v/ £ ■ Vg t + . . . , (4.22) 

which has exactly the structure of Eqns. (4.5) - (4.11). 



For notational simplicity we will often drop the index l s' denoting short-wavelength quantities, 
but it is important to keep in mind that the fields entering the effective stress-energy for long- 
wavelength perturbations are the short-wavelength fields v s and 4> s . In the remainder of the paper 

will always refer to [t^] 8 + [r^f 2 ps [t^] 3 . 
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When considering the effect on scales with fc< Awe will often replace the Fourier modes 
inside the small domains defined by the scale A by their ensemble averages, 

(Ma)(x). (4.23) 

The ensemble average is used to determine what the spatial average would be in a typical domain. 
The effect of random statistical fluctuations is captured by the variance, 

Var([r^] A ) = (Mi) - ([r^] A ) 2 , (4.24) 

i.e. the variance quantifies the expected statistical variation of the spatial average between differ- 
ent domains. This goes to zero as the domain size becomes large, but can be important in some 
applications, see §D.2. 

In the following we will assume this smoothing, [...]a, and ensemble averaging, (...), in the 
long- wavelength limit, but we will often simplify our notation to 

dvk) -> M -> ( 4 - 25 ) 

The meaning of this notation should be clear from the context. 

Finally, we remark that the spatial dependence of Eqn. (4.23) can arise from two different 
effects: 

1. Long-wavelength fluctuations (Section 5): Long- wavelength density fluctuations correlate 
small-scale fluctuations in spatially separated regions. 

2. Random statistical fluctuations (Appendix D): Stochastic fluctuations arise from quantities 
in a given realization of the universe being different from their ensemble averages. 



4.2 Truncation of the Boltzmann Hierarchy 

In the next section we will interpret the coarse-grained effective stress-energy tensor {t^ u ) as 
an effective fluid with density, pressure and anisotropic stress. This fluid description relies on a 
truncation of the hierarchy of moments of the distribution function for dark matter particles. We 
therefore digress briefly to argue that this truncation naturally arises for collisionless dark matter 
particles in an FRW universe. Our discussion is based on the classic treatment of the Boltzmann 
hierarchy by Ma and Bertschinger [56]. 

A collisionless gas of non-relativistic particles is described by the particle phase space density 
f(x, p, 77), where fd 3 xd 3 p is the number of particles in an infinitesimal phase space volume. Here, 
x is the comoving spatial coordinate of the particle and p = amx its conjugate momentum. 
Phase-space conservation leads to the collisionless Boltzmann equation (or Vlasov equation) [56] 

+ — ■ V/ - omVI ~ = . (4.26) 
orj am op 

Moments of the distribution function / define the particle density, p, the particle momentum flow, 
pv, as well as isotropic and anisotropic contributions to the pressure, p and a (see Appendix A). 
Inhomogeneities and anisotropies are described perturbatively 

f( X ,p,r ] ) = f (p)[l + S f (x,p,p,r l )] . (4.27) 
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Let us assume that we are in the quasi-linear regime, and expand the Fourier modes of the 
perturbation Sf in terms of Legendre polynomials P n , 

oo 

S f (k,p,p, r?) = ^H) n (2n + l)5 [ ;\k,p, V )P n (k ■ p) , (4.28) 

n=0 

where k and p are vectors of unit norm. The first few moments of the perturbation foSf then 
give the perturbed energy density, pressure, energy flux, and shear stress 

Sp = ma' 3 J 4Trp 2 dpf (p)5 [ j ] , (4.29) 

1 F 2 

5p = ma - 3 Wdp^h{p)8f , (4.30) 

O J (X Tfl 

(p + p)e = ma- 3 k f 47Tp 2 dp—f (p)5 [ f ] , (4.31) 
J am 1 

(p + p)a = -\ma- 3 f ^p 2 dp ^f (p)Sf . (4-32) 

o J (X Tfl 

We now argue that higher moments are systematically suppressed for scales larger than the 
non-linear scale, i.e. 

6f~ 3] < 6f, for A;>A: NL . (4.33) 

First, we note that the Boltzmann equation (4.26) implies a hierarchy of evolution equations for 
the higher moments [56] 



6^ = kv p 



n + 1 § [n+l] _ n [ W -1] 

2n + l / 2n + l / 



n > 2 , (4.34) 



where v p = = x is the particle's peculiar velocity which is to be distinguished from the mean 
peculiar velocity entering 9 = V • v. At long wavelengths, any time evolution is of order the 
Hubble time Ti' 1 , so that Eqn. (4.34) may be estimated to give 



n>2. (4.35) 



5 [ P ~ kv v U- x ^-S f [n + l]--^—Sf 
J p 2n+l /L 1 2n + l * 

If kvpli- 1 is much smaller than unity, this leads to a natural hierarchy of higher moments 

Sf ] ~ (kvpH' 1 ) 71 - 2 Sf ] . (4.36) 

This is similar to what happens in a conventional fluid, for which an equation like Eqn. (4.35) 
applies with prefactor kv p T c , where r c is the characteristic collision time. In a sense, the presence 
of such a hierarchy is really the definition of a fluid. In our case, v p t c is the mean free path and at 
distances larger than the mean free path, kv p r c C 1, we can truncate the hierarchy and describe 
the system by coarse-grained fluid variables. In the absence of such a hierarchy, our approach for 
deriving a long-distance effective fluid would be doomed. 

For the dark matter fluid is not obvious that kv^H.- 1 should be small since k/H 3> 1 inside 
the horizon. To see that the truncation of the Boltzmann hierarchy indeed occurs we need to 
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estimate the typical particle velocities. As a bound on the maximal particle velocities we take 
the velocity at the non-linear scale 

v p < A v (k NL ) , (4.37) 



where (see Appendix B) 



This implies that 



A;(A- NL )~As(^ L )|^~|^- (4-38) 

NL NL 



kVpH' 1 < — , (4.39) 
which is indeed small for scales larger than the non- linear scale, k <C &nl- The Boltzmann 

[21 

hierarchy therefore can be truncated beyond 5^ (the moment responsible for anisotropic stress) . 
This may be understood intuitively: during a Hubble time particles don't move more than a non- 
linear distance, not because they interact strongly like in a conventional fluid, but because they 
move slowly and haven't had time to move far. The hierarchy between the Hubble scale Ji~ x and 
the non-linear scale plays a fundamental role in this interpretation. Gravity induces a finite age 
and a finite horizon, which allows for the formation of a hierarchy among the multipoles. On the 
contrary, in the absence of gravity, the system could be infinitely old, and we would be left with 
no hierarchy and therefore no fluid. It is in this sense that we refer to the result of integrating 
out short- wavelength non-linearities in a FRW universe as a 'gravitational fluid', highlighting the 
fundamental importance of gravity. 

4.3 The Effective Fluid 

Let us now begin to discuss the physical content of the stress-energy tensor arising from second- 
order matter fluctuations at short scales [k S> A) as described by Eqns. (4.8) and (4.9). We 
argued in the previous subsection that at long wavelengths (fc < A) the effective stress-energy 
tensor can be put into the form of a fluid with density, pressure and anisotropic stress, 

Pes = MUfU? « -<r°o) , 3 PcS = (r^hf ~ (t\) , (E^eff ~ (f^) , (4.40) 

where f ^ is the traceless part of t ^ , and we have absorbed a small correction from bulk viscosity — 
cf. Eqn. (3.14) — into the definition of p e ff- We notice from the defining equations of T Mi , that our 
description of the effective fluid requires as the basic input the following correlation functions of 
the short-wavelength fluctuations, 

Ktj = \{{l + 5)v iVj ) (4.41) 
a,- = - ( MA ~ (Mil ( 4 42) 

In the last passage of Eqn. (4.42) we have assumed that (4>4',i) = 0- As we will discuss in 
§5.2, this quantity is in general non-zero and proportional to a derivative of long- wavelength 
fluctuations. This term is therefore irrelevant in the long- wavelength limit. For the moment we 
neglect this effect and concentrate only on the long- wavelength limit. With the above definitions, 
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the spatially-averaged components of the effective stress-energy tensor, Eqns. (4.8) and (4.9), can 
be written as 



(A) = -p(1 + 2k-5w), (4.43) 
(ry = p(2Kij + uSij - 2ujij) , (4.44) 

where k and u are the kinetic energy and the potential energy of a fluid element normalized with 
respect to the background energy density p, 

n = Ku = -{(1 + 5)v 2 ) and u = w« = -(5 0) < . (4.45) 

The spatial averages of the kinetic and potential energies are related by the cosmic energy equation 
(or Irvine-Layzer equation) [23] 

d -^l + n(2 K + u) = 0. (4.46) 

Finally, we note that the properties of the effective stress-energy tensor, Eqns. (4.43) and 
(4.44), are expressed in terms of the power spectra of small-scale density and velocity fluctuations, 
e.g. 

'dlnq^Afa), (4.47) 

where the integral is over modes with q > A. This information is available through the halo 
model [55] or iV-body simulations [33, 34, 35]. Alternatively, the cosmic energy equation (4.46) 
relates k and w given some initial conditions [32]. 

4.4 Renormalization of the Background 

From the above analysis it is straightforward to see that integrating out short-wavelength fluc- 
tuations leads to a renormalization of the background. We define the new background as the 
k <C A limit of the effective fluid, 

M = " ITa {t ° o) ' 3pcS = ITa {A) ' ( ^' )cff = ITa { ^ } ■ (448) 

Eqn. (4.48) describes the fluid on very large scales, where spatial fluctuations are suppressed by 
k 2 with q* the typical scale of non-linearities. In particular, on superhorizon scales these 
fluctuations are highly suppressed. We defer a treatment of fluctuations on subhorizon scales 
until §5. 

Density. We find that the effective energy density receives contributions from the kinetic and 
potential energies associated with small-scale fluctuations 

p cff = p(l + 2k - 5u) . (4.49) 
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Introducing the 'boosted mass density' 17 (see Appendix B) 

p m = l{v)e~^p, (4.50) 
= p(l + y-3^ = P +- P (k-Quj), (4.51) 

the effective density of the fluid is 

PcS = (iPm) + ^(Pm4>} , (4-52) 

= p m (l + K + u). (4.53) 

This shows that the background energy density is corrected precisely by the total kinetic and 
potential energies associated with non-linear small-scale structures, in perfect agreement with our 
analysis in §2. For virialized scales with 2k + cj = the cosmic energy equation (4.46) implies 
that the sum of kinetic and potential energies, k + u>, is constant, so that virialized scales only 
contribute a time-independent renormalization of the background density. 

Pressure. The effective pressure of the fluid is 

3p cS = p rn {2K + uj) , (4.54) 

and its equation of state is 

W eS = ^ = 1(2k + u). (4.55) 
Pes 3 

We see that for virialized scales the effective pressure vanishes. As intuitively expected, a universe 
filled with virialized objects acts like pressureless dust. (This agrees with the conclusion reached by 
Peebles in [44].) Non-virialized structures, however, do have a small effect on the long-wavelength 
universe, giving corrections to the background of order the velocity dispersion, 0(v 2 ). We will 
have more to say about this in §6.4. Furthermore, in Appendix D we will show in perturbation 
theory that 2k + uj > {e.g. in linear theory 2kl + ujl = \^~L > in Einstein-de Sitter). Peebles 
[23], using the cosmic energy equation (4.46), shows that the inequality 2k + uj > even extends 
to the non-perturbative regime. We have therefore shown that the induced effective pressure is 
always positive, p e s > 0. 

Anisotropic stress. On very large scales the anisotropic stress (£^) e ff averages to zero, i.e. it 
has no long-wavelength contribution: 

lim(£%) eff « . (4.56) 
fc<A J 

This straightforwardly follows from the isotropy of the fluctuation power spectrum. On very 
large scales, the gravitationally-induced fluid therefore acts like an isotropic fluid; its only effects 
are small 0(v 2 ) corrections to the background density and pressure. Anisotropic stress, however, 
does become important when studying the evolution of perturbations on subhorizon scales. We 
discuss this limit in §5. 



17 The physical role of p m becomes clear when thinking about dark matter in terms of particles and their conserved 
number density, n oc p m . A relativistic boost relates the physical number density in the inertial frame where the 
fluid velocity is v to the number in the fluid rest frame, n rest oc p. The factor j(v)e~ 3 ' 1 ' in Eqn. (4.50) clearly relates 
the volume element in the rest frame to the physical volume in the moving frame. 
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4.5 Decoupling and the Virial Theorem 



In the Einstein-de Sitter universe, VL m = 1, the Newtonian potential is scale-invariant on all 
scales and the main contribution to the effective gravitational stress-energy comes from very 
short-wavelength modes. 18 In principle one may worry that the effect becomes arbitrarily large 
for modes with k S> %. In this section we show that modes of very short wavelengths, in fact, 
give no effect if they correspond to virial equilibrium. 

In the previous section we have already seen an example for the decoupling of virial scales in 
the fact that they do not contribute an effective pressure. Here, we will demonstrate that this 
result extends to the full tensorial version of the virial theorem. Specifically, we will prove that 
the spatial part of the effective stress-energy tensor, averaged over a domain of size A -1 , is equal 
to the second derivative of the inertia tensor associated with that domain, 

M* = lw- (457) 



Virial scales decouple from the long-wavelength dynamics. 
Proof I: 
Consider 

[nM x ) = / W A (\x-x'\)T t:j (x') (4.58) 

J x' 

= [pViVj] A - [<)>,k<t>,Aj - 2(j) ti <j>j] K . (4.59) 

Since we are interested in the effect of virial scales deep inside the horizon we ignored the expansion of 
the universe by setting a = 1. There is a subtlety associated with this assumption that we will discuss at 
the end of the next proof. The first term in Eqn. (4.59) may be written in terms of the tensor 

K tJ = \J f W K {\x - x'\) pv iVj (x') , (4.60) 

whose trace, K = Ka, is the non-relativistic kinetic energy of the fluid. We now show that the second 
term can be written as the tensor 

Wij = \J W A (\x - x'Dpix'^ix') , (4.61) 

where 

$y(aO = -G [p(x' + r) T -^. (4.62) 

J r ' 

The trace of Wij is the gravitational potential energy of the fluid. 

18 In the Einstein-de Sitter universe the non-linear scale &NL is always a fixed fraction of the horizon scale, 
ktih/'H = const. In the real universe, the scale of matter-radiation equality is always bigger than the non-linear 
scale. The effect of non-linear scales therefore receives an additional suppression from the fact that they entered 
the horizon during radiation-domination and then decayed until matter-domination. In that case the biggest 
contribution is expected from modes with k ~ k eq [13]. 
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Consider the Fourier transform of the term (p^j, 

/ ' ' k " J 'o.,o.. l = - / q\k - q) j (t> q (j)k- q 

J as J a 



-+ (4ttG) 2 / ^- PqPk - q . 
J q y 



(4.63) 
(4.64) 



In the second line we used the Poisson equation in the Newtonian limit, q 2 (f> q = 4nGp q , and took the limit 



k -C q. For the gravitational term in Eqn. (4.59) this implies 

* = -g^j J e~ lkX {(t>,k(t>,k5ij -2(j) ti <j)j) = ~ J ^[q 2 5 lj -2q i q j \p qPk ^ 



Using the identity 



and hence 



we can write Eqn. (4.65) as 



iq .x<m _ _L 



tin - 



13 'xl 2 



/ e^[q%-2q^] = 
J a y 



,\j / ,1 I 



Pk- 



q ' 



(4.65) 

(4.66) 
(4.67) 

(4.68) 



where &ij(q) is the Fourier transform of Eqn. (4.62). Considering the Fourier transform of Eqn. (4.61), 



Wa(k) = \w A {k) ■ J ^M)Pk- 



Q 1 



wc find 



Hence, 



8ttG 



[4>,k4>,kkj - 20,i</>,j]A = Wij . 



] A = 2Kij + Wij . 

In a background with non-zero pressure p we would have found 



[T ij ] A = 2K ij +W ij +P5 ij , P= f W AP . 



(4.69) 

(4.70) 
(4.71) 

(4.72) 



It is a standard result from hydrodynamics {e.g. Ref. [57]) 19 that the r.h.s. of Eqn. (4.72) equals the second 
time-derivative of the moment of inertia tensor 



Jy = / W^X-x'Dpx'tX'j 
Jx' 

Hence, we have proven the desired result, Eqn. (4.57). 



(4.73) 
QED ■ 



9 In the next proof we will not actually need to rely on standard results from fluid dynamics. 
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Virial scales decouple from the long-wavelength dynamics. 
Proof II: 

Eqn. (4.57) also follows straightforwardly from conservation of the total stress-energy tensor t m „ for 
matter plus gravity, 

d^ v = , (4.74) 

as we now show. For the moment we are ignoring the expansion of universe, since we are interested in 
scales well inside the horizon. Consider then the following integral: 

0=[ d k d l [W A (\x-x'\)x i x j T kl ] . (4.75) 

J x' 

This vanishes because the integrand is a total derivative, and it is localized. Expanding the derivatives we 
get 

= f W A (\x - x'\) hr 13 + x { x j d k diT kl + 2 (x'Si + x j Si)d k T kl ] , (4.76) 

J x' 

where (and henceforth) we neglect derivatives acting on W\, because we are assuming that the smoothing 
scale A -1 is much larger than the typical variation scale of t^ v ' . Upon integration by parts, the last term 
can be written as minus twice the first term. The second term instead, can be converted into a second 
time-derivative by using conservation of t^ u twice. We finally get 







/ W A (\x-x'\) [-2r^+a 2 (xVr 00 )] , (4.77) 

Jx> 



which is exactly Eqn. (4.57). 

Notice that this proof holds at all orders in the post-Newtonian expansion, since Eqn. (4.74) does. 
So, even if we consider a relativistic virialized system, say a black hole-black hole binary, the total stress 
tensor r^, smoothed over a scale larger than the system's size, obeys Eqn. (4.57). 

Now, the expansion of the universe is negligible at subhorizon scales, but can its effect build up over 
time and give 0(1) modifications to Eqn. (4.57) on time-scales of order Hubble? To show that this cannot 
happen, it is particularly convenient to use the Newtonian coordinates of §2. The background FRW metric 
is given in Eqn. (2.1). Consider then a virialized system localized around x = 0, of size I much smaller 
than Hubble. In these coordinates the effect of the FRW expansion on this system is manifestly a tidal 
effect suppressed by H 2 £ 2 , with nothing explicitly becoming large 20 for t <~ i/ -1 . Of course, even such 
a tiny tidal-coupling does perturb the system to some extent, and it will have an integrated effect over 
long time-scales. However, for this effect to become of order one, one needs to wait for a time of order 
if _1 (iJ _1 /f s ) ^ if -1 , where t s is the system's dynamical time. Notice that this is perfectly compatible 
with there being the usual secular redshift effects associated with the Hubble expansion. For instance 
right at x = tidal effects are strictly zero, yet the background p does redshift according to the standard 
FRW equation 

p = -3H(p + p), (4.78) 

which yields an order-one variation of p over a Hubble time. In our coordinates it is clear that this is not 
a gravitational effect. Rather, it is enforced by the specic^-relativistic continuity equation, applied to the 
Hubble velocity field v ~ Hx. 

In conclusion, the expansion of the universe does not perturb Eqn. (4.57) appreciably, even on time- 
scales of order Hubble. QED ■ 



20 In fact for spherical systems there is no effect at all, due to Birkhoff's theorem. So the effect is really a 
tidal-coupling between the background expansion and the traceless quadrupole moment of the system. 
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Note that for virialized sources the r.h.s. of Eqn. (4.57) is not just zero. Rather, it is of the 
same order of magnitude as for non- virialized sources, namely ~ pu 2 . The crucial property of 
Eqn. (4.57), however, is the explicit time-derivative structure. Since we are interested in the 
effect of short-scale structures on the long-scale and long-time evolution of the universe, we are 
effectively only sensitive to the average of Eqn. (4.57) over a long time T, much longer than the 
typical short-scale dynamical time t s , 



d 2 Iij \ 1 dhj 7 



drf It T drj 



o T 



l fP vir, pv 2 (^) <^pv\ (4.79) 



where t is the typical size of virialized structures. Equivalently, in Fourier-space we are interested 
in the very low frequency (a; ~ H) component of Eqn. (4.57), which is automatically small because 
of explicit powers of the frequency, 

d 2 h 



. (4.80) 



drj 2 

Note that this Fourier-space argument does not go through for non-viralized sources, because their 
motion is peaked at low frequencies. For instance, for a free particle with p = m5T)(x — xo{t)) 
and xo(t) = Vot we have 

Iij(t) = mviVjt 2 5r,(x - x (t)) , (4.81) 

whose second time-derivative is unsuppressed, and in fact supported exclusively, at low frequen- 
cies. 

This decoupling (or non-renormalization) result proves that we shouldn't expect large effects 
from small-scale non-linearities on the background evolution. The gravitational non-linearities 
naturally cut off at very small scales with the main contribution coming from q* ~ /cnl (in the 
Einstein-De Sitter universe) or q+ ~ k eq (in a matter-radiation universe). 



5 Effective Theory of the Perturbed Fluid 

We have seen that on very large scales, k <C H, the effective anisotropic stress induced by second- 
order scalar fluctuations averages to zero. On superhorizon scales the effective fluid is therefore 
well-approximated by a perfect fluid with no dissipation. However, on scales comparable to 
or smaller than the horizon, k > H, we expect higher-derivative terms to become relevant. 
For example, we will see that anisotropic stress becomes important for the evolution of long- 
wavelength perturbations. The properties of this imperfect fluid are treated in this section. 

In §5.1 we introduce the effective theory of the fluid in the presence of long- wavelength 
perturbations. The theory is characterized by the UV-IR coupling of short and long modes. 
In §5.2 we illustrate the main effects for a universe filled with cold dark matter particles. We 
emphasize the importance of tidal forces from the long modes on the dynamics of the short 
modes. We describe, in §5.3, how the basic parameters of the effective theory are determined in 
perturbation theory, from matching to numerical simulation, or by fitting to observations. In §5.4 
we discuss corrections to the results arising from stochastic fluctuations of the fluid parameters. 
Technically these fluctuations arise because quantities in a given realization of the universe are 
different from their ensemble averages. Finally, in §5.5, we give a brief summary of our results. 
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5.1 The Effective Theory 



The basic degrees of freedom of the long-wavelength theory are the density contrast bi and the 
velocity divergence 21 Qi = V • V£. The (linearized) evolution equation for the velocity divergence 
is 

i 

(5.1) 



2 Pi 



In this section, we will argue that the non-linear UV-IR coupling of the perturbations leads to the 
short-wavelength source term (nj) being a function of Se and 0£. Intuitively, this arises because 
the presence of long-wavelength modes changes the particle geodesies on small-scales, leading to 
small changes in the stress-energy associated with short- wavelength modes. These effects originate 
from tidal forces and are hence proportional to derivatives of the long- wavelength gravitational 
potential (in §5.2 we will give an explicit example of how these effects arise in a dark matter 
universe). This leads us to introduce an effective approach, where we expand the source term in 
Eqn. (5.1) in derivatives of the long-wavelength fields, 



(nj) = p 



Cl 



(v 2 s )d 2 



+ p 



(4 



U 2 

n) 



C2 



*3 
CP 

A 2 



u 2 



+ 



+ 



(5.2) 



A 2 J 



+ •••){«!, 5e<f>i , ■ ■ ■ } 



Here, q and d^ are time-dependent coefficients of order one. The index 'n' on d^ runs over the 



fields represented in the vector between curly brackets, {v 2 , 



} . The double index jj in each 



term comes from the derivative-structure of the individual terms (constructed from combinations 
of didj, do and 5ijd 2 ). Similar expansions hold for Too and ro«. 

The terms in the first line of (5.2) are a derivative expansion in powers of k 2 /k^ L . It is easy 
to see how this scaling arises if we think about dark matter in terms of particles: the particles 
feel tidal forces from gradients of long-wavelength gravitational modes. Since the particles travel 



a distance of the order of v s H 1 



k' 1 



during a Hubble time, this leads to the k/k^L scaling. 



Alternatively, we may recall the description of dark matter particles via the Boltzmann equation 
in §4.2. There we saw that the moments of the distribution function are connected by a hierarchy 
of equations. Integration of the Boltzmann equation leads to a hierarchy of moments, with 
higher moments being suppressed by powers of fe/^NL- These higher moments feed back into the 
equations for the first two moments: the density and the velocity. The net effect is equations for 
the first two moments with additional higher-derivative terms. In an effective approach, we simply 
introduce these higher-derivative terms directly into the equations of motion and parameterize 
our ignorance of the theory for the short modes by the coefficients of these higher-derivative 
terms. In the next subsection we shall show explicitly how these terms arise in a universe filled 
with only dark matter. 



21 Below we will also comment on the vorticity »f = Vx ve. In fact, from a Lagrangian/Hamiltonian viewpoint 
the mechanical degrees of freedom parameterizing a fluid's configuration space are three scalar fields ^(Xjt), 
I — 1,2, 3 — the internal coordinates of the fluid element located at x at time t [21]. However, at the level of the 
equations of motion the density and velocity fields suffice. 
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The terms in the second line of (5.2) represent additional higher-derivative terms that are 
suppressed by powers of k/A, with A -1 being the smoothing scale. These terms arise from 
smoothing out the short modes to obtain equations for the long modes. Contrary to the terms 
arising from tidal forces, these terms need not be proportional only to the gravitational potential, 
but, at least in principle, can be proportional to any of the fluctuating variables. These fields 
are represented in the vector {vj,5e4>e, ■ ■ ■}■ We point out that the terms originating from the 
smoothing have to start at quadratic order in fluctuations - see e.g. Eqn. (4.10). Furthermore, 
since for a non-relativistic system 5g appears at non-linear level only when accompanied by other 
fluctuating fields, terms like <5|, <5|, ■ ■ ■ do not appear. This makes these higher-derivative terms 
quite subleading. 

Notice that all of the corrections to (r^) go to zero as k — > outside of the horizon (once 
the zero mode has been subtracted; see §4). It is convenient to use the linear Einstein equations 
to express 4>£ and its derivatives in terms of 5e and 9g. Concentrating only on the leading terms 
arising from two- derivative terms like didj(ftg, we obtain 

hkj (Tjj) 2 2 0£ 

-p y - c s d e - c vis — + ■ ■ ■ . (5.3) 

We see that the effect of the higher-derivative terms is to induce in the long-distance fluid an 
effective pressure perturbation, parametrized by c 2 , and an effective anisotropic pressure (or 
viscosity), parametrized by c 2 is (see §5.3). Notice that we parametrized the coefficient of the 
viscosity as c 2 is in dimensional analogy to c 2 , but the viscosity is of course a dissipative term 
that does not induce propagating sound waves. The following is worth noticing: c 2 and c 2 is are 
comparable in size and lead to comparable effects in the Euler equation, since, at linear level, 
OefH ~ de/T-L ~ Sg. This makes our effective fluid an unusual one. In conventional fluids the 
pressure is hierarchically larger than the viscous term by an amount of order k/k c , where fc" 1 
is the mean free path of the fluid particles. In our case the pressure happens to be very small, 
making it comparable to the viscosity. 22 

For simplicity, vorticity has been neglected in Eqns. (5.1)-(5.3). This was not strictly neces- 
sary. In practice, the fundamental description of our effective fluid is given by Eqn. (5.2), which 
describes the effective stress-tensor that enters in the Euler and continuity equations, without 
any assumptions about vorticity. In reality, since t^ v contains viscous terms and the fluid is 
compressible, even if we start with no vorticity at early times, it will be generated during the 
evolution (see Eqn. (5.28) below). This is an important difference to standard perturbation the- 
ory. In fact, in the standard approach, dark matter is described as a perfect and pressureless 
fluid, implying that vorticity cannot be generated unless the velocity field is defined through a 
proper smoothing similar to our definition (see for example [22]). Numerical simulations show 
that vorticity is generated when non-linear structures are formed. Our viscosity parameters are 
to be matched to simulations and we expect that the resulting fluid should be able to reproduce 
the generation of vorticity. 



22 This is due to the fact that the speed of sound is small, so that the time derivative of 5t is dominated by the 
Hubble flow rather than the effect of the gradient terms. Notice, however, that the effective stress-tensor is still 
that of a fluid. Higher moments of the particle distribution function are still suppressed with respect to the first 
two moments by factors of fc/fcNL- ft is this feature that ensures that the Boltzmann hierarchy can be truncated 
and that an effective fluid description is possible (see §4). 
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Ultimately, the parameters of the effective fluid in Eqn. (5.3) - or more precisely in Eqn. (5.2) 
- can be calibrated with numerical A-body simulations (§5.3), estimated in perturbation theory 
(Appendix D), or even fitted to observations. However, before discussing this, we would like to 
understand the effects more intuitively. 

5.2 Tidal Forces and Anisotropic Stress 

It is interesting to see how the higher-derivative terms of the effective theory arise in a concrete 
example. In this subsection, we consider a universe filled with cold dark matter only 23 and study 
the dynamics on very short scales. We will see explicitly that in this context, the higher-derivative 
terms arise because the long modes affect the motion of the short modes through tidal forces, 
and this then backreacts on the motion of the long modes themselves. To see this effect, we find 
it convenient to revert to thinking of the dark matter as a collection of particles. 

Consider the universe as a collection of small cubes of linear size A -1 . We are interested in 
the effects of long-wavelength modes on the dynamics inside each small cube. The equations of 
motion for each particle 'a' in the cube are 

lirf + = -V0s(*a) - V&(*«) , (5-4) 

where we have split the Newtonian potential into short and long modes, (fr s and <f>g, relative to 
the size of the cube. We decompose the particle positions into the center of mass of the cube, x, 
and coordinates relative to the center of mass, 5x a , with 



Y,Sx a = 0. (5.5) 



For simplicity, we have assumed that all particles have the same mass, m a = m, but the method 
generalizes straightforwardly to a distribution of masses. Summing Eqn. (5.4) over all particles 
in the cube gives the equation of motion for the center of mass 



d 2 x _^_^dx 
drj 2 dr\ 



2+ n — = - V Mx), (5-6) 



where we have used the fact that the sum over the short-wavelength modes averages to zero, 

£v&(aj a )«0. (5.7) 

a 

Eqn. (5.7) essentially represents the fact the internal forces taken into account in cf> s do not 
accelerate the center of mass. Subtracting Eqn. (5.6) from Eqn. (5.4) then gives an equation for 
the particle motion relative to the center of mass, 

^ = -V* - + -. , (5.8) 

drj z drj oxiOXj 



23 0ur approach of adding higher-derivative terms as in Eqn. (5.2) is of course more general than that: a universe 
filled with dark matter represents only a particular UV completion of our effective fluid. 
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where we have expanded the long- wavelength potential in a power series around x. We see that 
the effect of the long-wavelength mode <pi on the internal dynamics of particles in the cube is 
captured by the tidal tensor, didj4>£. 

Even in the absence of the long-wavelength potential the geodesic equation (5.8) typically 
cannot be solved exactly. In the Zel'dovich approximation Vi<p s ( x a) is evaluated at some initial 
time 770, when the particles have initial positions Xq [23]. The lowest-order approximation to the 
particle deflections then is 

5x l a ( V )^-b( V )V t MO, where %) ~ U~ 2 . (5.9) 
The tidal force term in Eqn. (5.8) is then included perturbatively by considering the source term 



' dx l dxi 



-m-Z^hW.- ( 5 - 10 ) 



Eqn. (5.10) reveals explicitly the UV-IR coupling between long- wavelength modes and short- 
wavelength modes. In principle, the Newtonian potential can be a function of time, so the 
solution of Eqn. (5.8) with the source term (5.10) can involve an integral over time. However, 
during the matter era <fi is constant and the solution of Eqn. (5.8) takes the form 

d 2 <f>£ 

Sx l a (rj) « -b{ri)di<j) s - g{rj) q iQ^ i^* > where 9(v) ~ • ( 5 - u ) 
The associated three- velocities of the particles then are 



^ WCl ^~ +C2 — %3 — ' where c i,c 2 ~0(l). (5.12) 

For purposes of illustration, we now consider the term pviVj C r^- and average over short- 
wavelength fluctuations, 

(pv iVj ) ~ S^ + ^pipv!), (5.13) 

~ Mp^> + ^(p"2>. ( 5 - 14 ) 

where we used 

(di4> s dj(j) s ) ~T-L 2 5ij{v 2 s ) and didjfa ~ Hv*^ . (5.15) 

Here, we have gone back to the continuum description of the dark matter and dropped the particle 
indices a. The first term in Eqn. (5.13) describes isotropic pressure and arises even in the absence 
of the long-wavelength mode, while the second term gives anisotropic stress from the breaking of 
isotropy induced by the long-wavelength mode. It is easy to estimate that the next term from 
tidal forces scales as (Vg)d fie/H , in agreement with what we anticipated in Eqn. (5.2). As 
expected, the tidal effects are suppressed on superhorizon scales, where the Newtonian potential 
4>e is constant. However, the above analysis shows that on subhorizon scales tidal forces induce 
an effective anisotropic pressure of the form 

(EyJdf-^oJ-. (5-16) 
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where is the shear tensor defined as 

°i3 = V (i,j) ~ V k,k ■ (5-17) 

In addition, the trace of Eqn. (5.13) leads to an isotropic pressure perturbation proportional to 
5e and to a bulk viscosity proportional to 0£. Our effective theory, by construction, describes all 
these effects, as well as higher-derivative corrections. 

It should be emphasized that the tidal effects we described in this section are automatically 
captured in the results of iV-body simulations. Those numerical works solve the particle geodesic 
equations exactly (and in the presence of long- wavelength perturbations) and hence include tidal 
effects and anisotropic pressure. It should also be emphasized that this is not a general relativistic 
effect, but is completely captured by the Newtonian dynamics if the Newtonian equations are 
solved exactly. Nothing is therefore lost when the parameters of the effective theory are matched 
to the results of iV-body simulations. 



5.3 UV Matching 

We now describe the UV matching of the coefficients of the effective theory to either numerical 
simulations or directly to observations. For purposes of illustration, let us return to Eqn. (5.3), 
the expansion of the source term in the Euler equation to leading order in 5c and Qe = ^ (the 
dimensionless velocity dispersion), 

y-^-4A. (5.18) 

To extract the speeds of sound from Eqn. (5.18) we correlate A with long- wavelength perturba- 
tions 5i and G>£, 

A 5 = (5 e A) = c 2 s P S5 - c 2 vis P S e , (5.19) 
A e = (Q t A) = c 2 s P se - c 2 vis P ee . (5.20) 

In linear theory, 5^ = — and Pss = Pee = \Pse\i implying that this approach would only 
allow us to obtain a linear combination of the sound speeds, c 2 + c^ is . However, non- linear 
corrections will break this degeneracy, Pss / Pee / I Pse \ , and allow us to measure (or calculate) 
c s and c v i s separately, 



2 Pee A s - Pse A e _ , 

° s= PssPee-Pl ' (5 - 21) 



2 _ Pse As - Pss A e 

Cvis " PssPee-Pl ' { ' 

The parameter c v j s may be related to the standard phenomenological parameterization of 
imperfect fluids whose stress tensor may be written as 

Tij = puiUj + {p- (6)~fij + Sjj , (5.23) 
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The bulk viscosity ( and account for viscous terms. At leading order in an expansion in 
spatial derivatives, an ansatz for S^- that is consistent with the second law of thermodynamics 
is [18] 

Eij- = -r)Gij . (5.24) 

Here 77 > and C > are the coefficients of shear and bulk viscosity, respectively and Oij was 
defined in Eqn. (5.17). With Eqn. (5.24) the Euler equation becomes the non-relativistic Navier- 
Stokes equation. From Eqn. (5.24) we see that the viscous speed of sound is related to a linear 
combination of the shear viscosity 77 and the bulk viscosity £, 



To extract the shear viscosity of the fluid (and hence via Eqn. (5.25) also the bulk viscosity) 
we consider 

k z p 3 p 

Correlating B with 6^ gives the shear viscosity 77, 

B e = {Q i B) = -\ri^P e e. (5.27) 
P 

We remark in passing that viscosity is a source term in the generation of vorticity. We see this 
from the linearized equation of motion for we = V x ve, which, at leading order in derivatives, is 

w\ + Uw\ = \^V 2 w l £ (5.28) 
2 pt 



+ 0(8 2 „ »?,...). 



The terms in the second line represent sources for vorticity. Notice that they vanish if the fluid 
is incompressible. 

Finally, we may also consider correlations with the trace of , 

c = ^ = c ^-c?e,, (5.29) 

3p p 

which contains information on isotropic pressure and bulk viscosity. 

One-loop perturbation theory. As we observed above, a linear (or tree-level) calculation in gen- 
eral does not allow an extraction of all fluid parameters, unless it is assumed that the contribution 
from the bulk viscosity is parametrically smaller (see Appendix D) 

C| « c 2 s , . (5.30) 



37 



If we don't want to make the assumption (5.30) we have to go beyond linear order 
theory to compute c s and c v - ls separately via Eqns. (5.21) and (5. 22), 24 where 



in perturbation 



Pes = PL(k) + P^ 1 + , (5-31) 
Pee = PL(k) + P$ 2) +pj£' 3) , (5.32) 
Pse = PL(k) + P£' 2) +P& 3) . (5.33) 

Here, Pl(&) is the linear power spectrum and P^') are higher-order (one- loop) corrections. Ex- 
plicit expressions for the one-loop corrections p(* J ) may be found e.g. in Appendix A of Ref. [20]. 
From those expressions one may infer that the velocity fields are more sensitive to tidal fields than 
the density fields [41]. In fact, the non-linear growth of 9 is smaller than the linear prediction, 
i.e. on the relevant scales Peo(k) is suppressed relative to the linear P(k). The one-loop calcula- 
tion is feasible, but presenting it here would not be very illuminating. Furthermore, as Ref. [20] 
recently demonstrated, on small scales one-loop perturbation theory becomes increasingly unre- 
liable in capturing even the correct qualitative features of the velocity-density cross-correlation. 
Instead, we now describe how the parameters of the effective fluid should be determined from 
numerical simulations or directly from observations. 

N-body simulations. An A-body simulation of moderate box size allows one to measure the 
statistics of small-scale fluctuations and construct the stress tensor r^-. Furthermore, at the 
matching scale A the following correlations can be measured: P$s, Pse, Pee- Measuring the 
correlations of kikjTij with 5^ and 0g then provides a measurement of c s and c v i s via Eqns. (5.21) 
and (5.22). The basic idea is that a relatively small simulation suffices to normalize the fluid 
parameters and this fluid is then used in a perturbative treatment of long- wavelength fluctuations. 

Fitting to observations. Alternatively, the fluid parameters may simply be retained as free 
parameters of the effective theory. Computations are then performed with these arbitrary coef- 
ficients. Finally, the parameters are determined by comparing the results for certain observables 
to the corresponding measurements of the matter fluctuations. Notice that this may allow one to 
take into account effects not normally included in simulations (such as the effects arising from the 
presence of baryons), by fitting observations to a general one-component (or multi-component) 
effective fluid for the total matter density with free parameters as in Eqn. (5. 2). 25 

5.4 Corrections from Statistical Fluctuations 

We have used ensemble averaging and the ergodic theorem to estimate average quantities in 
domains of size A -1 . For any specific realization of the universe there will be a random statistical 
error in that estimate. In Appendix D we estimate the size of that effect in perturbation theory. 
Here, we discuss the basic conclusions of that analysis. 

24 0ne may ask how non-linear perturbation theory knows about dissipative effects. Even if the initial conditions 
are purely in the growing mode, non-linear interactions create both growing and decaying modes [24] . The decaying 
modes running in loops of the perturbation expansion carry the information on dissipative effects. 

25 This proposal is limited to sufficiently large scales for which relativistic effects, such as those due to the epoch 
of reionization, are small. 
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Consider for instance the effective pressure 



Pes = Pes + ¥eS + <Wff ■ (5-34) 

Here, p e s ~ c 2 s p denotes the renormalized background pressure (§4.4), 5p e eS ~ c%p8e is the per- 
turbation induced by long- wavelength modes Sg (§5.2), and Sp^^ is the perturbation induced 
by random statistical fluctuations, which enters in the equations for the long- wavelength mode 
as an external source. The statistical properties of this random source could be computed in 
perturbation theory (when this is a good approximation), or measured from simulations. In [42] 
the treatment and inclusion of these kind of terms has been developed in a different, but related 
context. In Appendix D we establish the following two facts: 

1. Both 5p e cS and Spf^ are much smaller than p e g if the momentum integral defining p e g is 
UV-dominated (i.e. if the integral is dominated by scales with q± » A, where, as before, 
g* is the typical scale of non-linearities: q* ~ ^nl in Einstein-de Sitter or q* ~ /c C q in 
our universe). Quantitatively, the stochastic fluctuations and the background pressure are 
related as follows 

^cff at = aPcfr, (5.35) 
where a is a random variable with variance 




(5.36) 



For A < q± stochastic fluctuations are therefore suppressed. This result may be understood 
as the standard ^= suppression of random fluctuations, where N is the number of space 
domains sampled. 

2. At the scale dominating gravitational non-linearities, q+, the two contributions to pressure 
fluctuations, Sp f cS and Sp^ 1 , are of the same size. On larger scales, (5p^ t , in fact, dominates 
over 5p £ eS — i.e. dp^ > 5p f cS — with equality holding at q+. However, at least at leading 
order (see [42]), the fluctuations 5pfg t are uncorrelated with the long- wavelength density 
fluctuations Sg, so their effect on the power spectrum of density fluctuations is subdominant 
for scales with k <C q*. On all scales, AP$ > Ai-*| tat . Nevertheless, for specific applications 
of the effective theory (see §6) the corrections from stochastic noise can be significant and 
will have to be included. 

We emphasize that stochastic fluctuations in our fluid are much larger than the corresponding 
effect in conventional fluids. In practice, the correlations of stochastic noise contributions there- 
fore have to be measured in the simulations or must be included phenomenologically as external 
sources. 

5.5 Summary 

Let us briefly summarize our conclusions so far: 

We have argued that the result of integrating out small-scale non-linear fluctuations can be 
represented at long wavelengths by an effective fluid. Studying the Boltzmann hierarchy for 
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dark matter particles, we have shown that the fluid is somewhat unconventional. The hierarchy 
is truncated at the second moment not because there are sizable interactions, but because the 
universe has finite age and particles only move a small distance (of order the non-linear scale) 
over a Hubble time. Higher moments therefore haven't had enough time to develop rather than 
being suppressed by interactions. Also unlike a conventional fluid, our effective fluid can have 
large stochastic fluctuations. 

In §4 we studied the renormalization of the background density and pressure induced by 
kinetic and potential contributions from small-scale structures. We proved that the contribution 
to the effective pressure is always positive and vanishes for viralized structures. 

In §5 we then studied the dynamics of long-wavelength fluctuations of the effective fluid. We 
have shown that dissipative effects, such as an effective viscosity, are crucial to understanding 
the evolution of long-wavelength modes. We furthermore introduced an effective theory for the 
perturbed fluid. Its key elements are: 

Hierarchy of scales. In the simplest application of the effective theory we have in mind the 
following hierarchy of scales (see Fig. 2) 

q* > A > k , (5.37) 

i.e. the characteristic scale of non-linearities g* is much smaller that the smoothing scale A, which 
in turn is smaller than the scales k to which we apply the theory. It is the large hierarchy between 
g* and k that allows for a perturbative treatment and a fluid description. The hierarchy between 
k and A ensures that higher-derivative corrections (see Appendix A) are negligible. The hierarchy 
between g* and A suppresses the effect of stochastic fluctuations (see Appendix D). 



■ 



conventional perturbation theory 
1 1 1 ► k 

effective theory A ft ^ 



Figure 2: Hierarchy of scales and perturbative expansions. The effective theory only contains modes with 
k < A, while conventional perturbation theory contains modes with k ~ /cnl where the perturbative 
expansion is known to break down. 



UV matching. Above the momentum scale A, and especially near the scale of non-linearities 
g*, the perturbative treatment of fluctuations breaks down. We therefore propose to numerically 
compute the effect of these non- linear scales. The cummulative effect of these small scales defines 
the parameters of the effective fluid (w e s, c s , c v i s , etc.) at the quasi-linear scale A: e.g. the 
parameters c s and c v j s appear as time-dependent expansion parameters in the source term for 
the Euler equation 

= *.h-&% (5-38) 

The more precise version of the effective theory given in Eqn. (5.2) contains further fluid pa- 
rameters. In §5.3 we described how these parameters can be extracted in A^-body simulations. 
Alternatively, the fluid parameters may simply be retained as free parameters to be measured by 
fitting perturbation theory calculations in the effective theory to the data. 
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Controlled expansion. After this matching calculation, the theory below the momentum scale A 
has well-defined small perturbation variables 5e(k <C A) and 6g(k <C A). On scales with k <C A 
we therefore define the theory as an expansion in 5i and 0g. The perturbation theory in Sg and 0g, 
evolving in a background characterized by w e s, c s and c v j s , is well-defined and under control. This 
is in contrast to other perturbative schemes, where the definition of small expansion parameters 
isn't as clear (see Ref. [20] for a discussion). 

RG flow. It is expected that the parameters of the effective theory like the speeds of sound 
c s (A) and c v i s (A) depend on the cut-off scale A at which the matching to numerical simulations 
is performed. In effective field theories this dependence of parameters on the matching scale 
is described by the renormalization group (RG) equations [5, 6, 7]. In our case, RG flow type 
equations will describe how the fluid quantities change as the cut-off is changed from A to A'. 
The corresponding beta-functions can be computed perturbatively in the effective theory. 

Loops and higher-derivatives. In principle, perturbation theory in our effective theory allows for 
arbitrary precision by going to higher and higher order in bt and 0£. This is different from the 
way in which most conventional cosmological perturbation theories are organized. In practice, at 
least the following two complications have to be considered: 

i) Going to higher order in perturbation theory requires also that the matching calculation in 
Eqn. (5.38) is performed to higher order in 5e and 0£. This is not a problem in principle, 
but it complicates the effective theory by adding extra fluid parameters — see Eqn. (5.2). 

ii) At some order, the higher-derivative corrections from the smoothing may not be negligible 
anymore and should hence be included in the effective t^ u — cf. the second line in Eqn. (5.2). 
At what point exactly these higher-derivative corrections become important will depend on 
the specific application and the hierarchy of scales it involves. 

6 Applications of the Effective Theory 

In this section we suggest possible applications of the general formalism we introduced in this 
paper. Our discussion will be rather schematic, leaving further calculations and numerical inves- 
tigations to future work [60]. 

6.1 An Alternative to Conventional Perturbation Theory 

In typical applications of effective theories one is interested in calculating the effects of some 
short-distance (high-energy) physics, characterized by a scale q+, on the dynamics at a long- 
distance (low-energy) scale k <C (/*. If the full theory is known, then these effects can in principle 
be calculated explicitly. If the full theory is not known (or the full theory is known, but not 
computable, e.g. strongly coupled) then high-energy effects are parameterized as k/q* corrections. 
Allowing all corrections consistent with a given set of symmetries accounts for the UV physics in 
a completely model-independent way. The parameters of the effective theory are then constrained 
either from numerical simulations or measured by experiments. 
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In this work we have applied the effective field theory philosophy to cosmological perturba- 
tions. The full theory (General Relativity coupled to matter) is known, so in principle all effects of 
short-scale non-linearities can be computed explicitly. The parameters of our effective theory are 
determined either from iV-body simulations, or matched to observations. After the UV match- 
ing, the effective theory at long-wavelengths is that of an imperfect fluid characterized by only a 
few physical parameters like the equation of state w e f[, the speed of sound c s , and the viscosity 
c v is. On large scales the effective fluid has small perturbation variables (e.g. 5e < 1), allowing 
for a controlled perturbative expansion at low momenta. Conceptually, our approach therefore 
offers a well-defined and controlled treatment of the effects short-distance non-linearities on the 
long- wavelength universe. This is to be contrasted with the breakdown of most conventional 
cosmological perturbations theory techniques at small scales (see Carlson et al. [20]). Often the 
conventional approaches are formulated as an expansion in 5, with loop integrals running over 
modes with S S> 1. Furthermore, the equations that are solved are the equations of a pressureless 
fluid (but see [43]). Since both the equations and the perturbative expansion is known to break 
down at the non-linear scale, errors are expected. However, the size of those errors is often hard 
to estimate a priori. In contrast, at low momenta our effective theory has well-defined small 
expansion parameters, which solve equations that are corrected by the properties of the effective 
fluid. Furthermore, like in all effective theories, the errors can be estimated and are under control. 

Our examples for applications of the effective theory will work their way from inside the 
horizon to superhorizon scales: In §6.2 and 6.3 we describe how the fluctuations measured in 
baryon acoustic oscillations (BAO) and the perturbations arising in models of preheating are 
affected by the viscosity of the effective fluid. Then, in §6.4 and §6.5, we explain how the 
superhorizon evolution of curvature perturbations and the backreaction of non-linearities on the 
background FRW evolution are understood in the effective approach. 

6.2 Baryon Acoustic Oscillations 

BAOs as a probe of dark energy. Oscillations in the baryon-photon fluid, frozen into the 
matter distribution at decoupling, provide a standard ruler to constrain the expansion history 
of the universe (for a review see [61]). These sound waves imprint an almost harmonic series of 
peaks in the matter power spectrum Ps(k), corresponding to a localized feature in the two-point 
correlation function £$(r) at about 100 h~ 1 Mpc [62]. The position and shape of the BAO peak 
are powerful probes of dark energy. While the evolution of density perturbations in the early 
universe is linear and therefore easily calculated, at late times the non-linear evolution of matter 
makes the interpretation of the BAO peak more challenging. The non-linear clustering erases the 
oscillations on small scales, which shifts and broadens the peak [63, 64]. An accurate modeling of 
the effects of non-linearities on the BAO feature is essential in light of the high precision aimed 
for in future BAO surveys; e.g. a shift in the BAO scale of one percent generates systematics in 
the deduced dark energy equation of state parameter of about five percent [65] , which is not 
negligible compared to the expected statistical errors in the next generation of galaxy surveys. 

Damping by viscosity. It is suggestive to explain the non-linear broadening of the BAO peak 
by the viscous damping of density fluctuations in an imperfect fluid (see §5). In Fourier space, 
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the broadening of the BAO peak is described by the following phenomenological ansatz [63, 64] 

Pobs(k) = e-^ 2 P L (k) + P mc (k) , (6.1) 

where Ph(k) is the linear power spectrum and P mc {k) is a mode-coupling term. 26 The exponential 
damping of the linear power spectrum in (6.1) is precisely the effect we expect from the viscosity of 
the effective fluid, with the parameter E 2 being related to a time integral of a viscosity parameter. 
For instance for a viscous fluid in the absence of gravity with bulk viscosity £, sound waves get 
damped as ~ e~^ fe *. Related observations have been made in perturbation theory in Ref. [64], 
where the exponential factor in (6.1) arises from a renormalization of the propagator for the 
evolution of density fluctuations. Thanks to the effective theory approach, we conclude that 
measuring the broadening of the peak corresponds to measuring the viscosity of the effective 
long-distance fluid. We find this to be quite a powerful statement. 

Real-world complications. Although we defer a detailed analysis to future work [60], we here 
point out a few of the expected challenges in this specific application of the effective theory. 
First, we notice that the BAO scale /cbao lies between the non- linear scale /cnl an d the scale of 
matter-radiation equality k eq , 

k cq < k B AO < &nl • (6.2) 

Since the dominant contribution to the gravitational fluid comes from scales with q+ ~ /c cq , 
the hierarchy of scales described in §5.5 - see Eqn. (5.37) - cannot be exploited in the BAO 
application. Specifically, to study fluctuations at the BAO scale one has to choose the matching 
scale in between the BAO scale and the non- linear scale: /cbao < A < &nl- This implies that 
the integrals defining the fluid parameters at the scale A are IR-dominated (i.e. dominated by 
modes with k ~ A), and therefore the stochastic pressure discussed in Appendix D is large and 
has to be included. Furthermore, the non-linear scale and the BAO scale are not very different, 
and the expansion parameter in the long-wavelength effective theory is not very small (although 
it is still smaller than one, making the expansion well-defined). These complications will make 
the application of the effective theory to the BAO example a bit non-trivial. 

6.3 Preheating 

For dark matter perturbations in the late universe, the large hierarchy between the scale of non- 
linearities and the horizon scale makes the effects we discussed in this paper small. It would 
therefore be interesting to study effects in the early universe where this hierarchy is larger and 
the UV-IR coupling is potentially more significant. It is well known that at the end of inflation a 
dynamical instability called preheating can convert the energy of the homogeneous inflaton field 
into the excitations of other fields [95]. This process is highly inhomogeneous and non-linear and 
therefore mostly studied numerically in lattice simulations [96]-[101]. The scale of non-linearities 
during preheating is found to be as large as 0.01H -1 [99], compared to 10 H -1 , the scale of 
dark matter non-linearities today. An effective description of smoothed perturbation fields during 
preheating may therefore be of interest. 

26 The mode-coupling term P m c(fe) relates to the stochastic small-scale noise of §5.4. 



43 



6.4 Non-Linear Superhorizon Evolution of ( 

Inflation [66, 67, 68] famously produces a nearly scale-invariant spectrum of scalar perturbations 
[69]. These perturbations are commonly described by the perturbation to the spatial three- 
curvature £, where gij = a 2 (n) exp[C(r?, a;)] dij. To relate observations of the cosmic microwave 
background and the large-scale structure to the primordial value of any evolution of £ after it is 
created during inflation has to be taken into account [53, 54]. One can prove [70], that £ remains 
constant outside of the horizon if: i) all modes that contribute to the non-linear definition of £ are 
outside the horizon, and ii) the pressure is a unique function of the energy density p(p) (adiabatic 
condition). Each Hubble patch then evolves like a separate FRW universe and ( doesn't evolve 
in time. However, short-wavelength fluctuations with wavelengths much smaller than the Hubble 
scale evade this theorem and can in principle allow a non- linear evolution of £. This was recently 
emphasized in Ref. [2]. In this section we show that this non- linear superhorizon evolution of 
curvature perturbations is completely consistent with the 'separate universe' picture. We explain 
that short-wavelength fluctuations induce a non-adiabatic pressure in the effective fluid describing 
the long- wavelength dynamics. 

Superhorizon evolution in comoving gauge. Ref. [2] derived the second-order curvature pertur- 
bation in comoving gauge (T* = 0) as a function of its initial condition (o, 

C = Co - l^d-'didjiditodjto) ■ (6.3) 
This implies a time evolution of ( of the form 

C = -^5- 2 ^(^Co^Co), (6.4) 
where the source term on the r.h.s. can be written in Fourier space as 

(k ■ q)(k 2 - k ■ q) 



Co(g)Co(* - g) ■ (6.5) 



k 2 

Crucially, Eqn. (6.5) does not vanish in the k — > limit, 

«max „ A? (a) 

lim<S = -/ dqdQ q (k-q) 2 (6.6) 
k^o J ,v 2it 2 q 

and is dominated by the UV cutoff at large q. Notice that in principle we could imagine a universe 
in which density fluctuations are never large, the comoving gauge is well-defined on all scales, 
and the integral above can be taken to q max = +oo. There is then no way out: £ seems to be 
evolving. We will now argue that the apparent evolution of ( is due to the renormalization of the 
background density and pressure induced by the short-scale non-linearities. Once the background 
has been redefined to take this into account, the curvature perturbation ( defined with respect 
to this background is indeed constant on superhorizon scales. 

Superhorizon evolution in Poisson gauge. We first reformulate the problem in the context of 
our effective theory in Poisson gauge. In this way we will see that the apparent superhorizon 
evolution of £ has a simple physical interpretation. 
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Smoothing the i-i Einstein equation (3.23) we find 



U 
P 



T k 



+ 



J A 



where 



H 

P + P 



Sqe 



with 



diSq. 



(6.7) 



(6i 



The quantity Ci is the standard linear definition of the comoving curvature perturbation 27 [53]. 
Using the 0-i Einstein equation to replace 5qe, it becomes 



= <fr e — 



H^e + H^i) md 

n-n 2 



2<S> e + HV e 

3 n 



(6.9) 



where the limit applies during matter-domination (MD). On the r.h.s. of Eqn. (6.7) we identify 
the pressure and scalar anisotropic stress of the effective fluid, 



Pes = 



T k 



such that 



and 



po- c s = 



v- 2 didj ( t 



-13 



-5V 



u 



(PcS + PO-cs) ■ 



(6.10) 



(6.11) 



This is the standard evolution equation for the long-wavelength curvature perturbation Q for a 
fluid with non-adiabatic pressure p e g and anisotropic stress a c g [53] . Both the pressure and the 
anisotropic stress can be split into a zero mode (§4), perturbations induced by long- wavelength 
perturbations (§5) and stochastic terms (§5 and §D.2), 



Pes = p eS + Spls + 5pfi l , 



(6.12) 
(6.13) 



All terms in Eqn. (6.13) average to zero on superhorizon scales. Similarly, 5p f cS and dp^^ in 
Eqn. (6.12) go to zero on superhorizon scales. However, the non-adiabatic zero mode pressure 
p e s remains and sources a superhorizon time evolution of the curvature perturbation, 



lim Cf 

k^u s 



U 

P' 



rPctt 



(6.14) 



However, Eqn. (6.14) simply reflects the fact that short-wavelength non-linearities renormalize 
the background and therefore slightly change the FRW history of separated Hubble patches. 
The curvature perturbation ( defined with respect to this new background would indeed be 
constant on superhorizon scales. The superhorizon evolution of ( found in [2] has therefore been 
demystified. 

Comments on the effective spacetime. As we have just seen, not all gauges are created equal. 
To illustrate the differences in gauges for the purposes of defining an effective theory on large 

27 In fact, the comoving curvature perturbation is often denoted TZ, reserving ( for the curvature perturbation 
on uniform density hypersurfaces. On superhorizon scales TZ and £ are equal (up to a sign). 
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scales, let us consider the long-wavelength limit of the effective spacetime. We split the covariant 
components of the metric into long- wavelength fields, g^, and short-wavelength fluctuations, g^ u , 

9^=9^ + 9%- (6.15) 

We then notice that smoothing the inverse metric g^ v is not the same as the inverse of the 
smoothed metric g^ u , with corrections being of order (see Appendix B) 

9ap[gr9f] A - (6-16) 

Importantly, Eqn. (6.16) in general doesn't vanish in the limit k <C A. However, as we discuss in 
detail in §B.3, in a gauge where 

ITa 9 ^ y> £f A ^r^] A (6-17) 

the smoothed spacetime is well-defined up to negligible corrections. In that case the long- 
wavelength metric, g^ v , is well-defined and can be used to construct covariant actions such as the 
action of a scalar field if living in this effective spacetime 

V~9e gT 'd^ipd^ip . (6.18) 

As we show in §B.3, the limit (6.17) is satisfied in Poisson gauge, but is violated in comoving 
gauge. In comoving gauge the notion of an effective spacetime after integrating out small scales is 
therefore ill-defined. This explains why the backreaction of small-scale fluctuations on £ is easier 
to interpret in Poisson gauge (see Appendix B for further details). 



6.5 (Non-) Acceleration from Backreaction 

It has been suggested (see e.g. Refs. [71]-[89]) that non-linear backreaction of short-wavelength 
fluctuations on the background evolution could explain the observed acceleration of the universe 
without the need for dark energy. In our work we see no indication of such a large effect. We 
have shown that small-scale inhomogeneities simply lead to a small renormalization of the FRW 
background. 28 In fact, we were able to prove that very small scales (in virialized equilibrium) 
decouple completely from the background dynamics (see §4) . This means that the effects of small 
scales are naturally cut off even in a pure matter universe, in which the effects of non-linearities 
can naively be expected to be large. As we stressed in the Introduction, this decoupling of 
small scales is more properly a non-renormalization theorem, and goes beyond standard effective 
field theory decoupling [3, 5, 6, 7]. Furthermore, we found that the induced effective pressure 
is positive, p e ff > 0, with a magnitude set by the velocity dispersion of small-scale fluctuations, 
w e s ~ 0(v 2 ) <C 1. We therefore do not see how the backreaction of small-scale non-linearities can 
possibly explain the acceleration of the universe. To arrive at this conclusion we had to assume 
that the spacetime can be described by the Poisson metric with small perturbations or equivalently 

28 To be clear, cosmological observations do not observe the background directly. Besides renormalization of the 
background, inhomogeneities in the matter affect light propagation (e.g. via lensing) and hence also the measured 
cosmological observables. For a discussion of these distinct effects see [90, 91, 92, 93]. 
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that small scales are well-described by the Newtonian approximation. This assumption has proven 
successful in many cosmological applications, but see Ref. [94] for a discussion of this point. 

Although we did not prove this, we believe that this decoupling of short-wavelength non- 
linearities should even apply in the extreme case that the universe is filled with a gas of black 
holes. In this case, our perturbative scheme breaks down. However, there is a large separation 
of scales between the horizon scale and the scale at which metric perturbation become of order 
one. 29 This means that one can match our effective theory to the effective theory of black holes 
by Goldberger and Rothstein [8]. In [8] it was shown that at intermediate scales, between the 
scale at which curvatures are of order one and the horizon scale, a collection of black holes 
may be described as point particles with some higher-derivative couplings to the metric. These 
particles can be taken to be the dark matter particles that we deal with in this paper, and 
our conclusions are therefore expected to apply straightforwardly even in this extreme example. 
Moreover, we proved that decoupling does hold for virialized systems containing black holes. We 
believe this argument makes a large backreaction from gravitational non-linearities even in this 
case impossible. 

7 Conclusions 

Cosmology has made significant progress by studying linear perturbations around a homogeneous 
Friedmann- Robertson- Walker background [102]. A crucial simplification of the linear treatment 
is that large scales (IR) are decoupled from small scales (UV), i.e. at linear order different Fourier 
modes evolve independently. Furthermore, in linear theory there exists a useful classification of 
the perturbations into independent scalar, vector and tensor modes. Beyond linear theory, the 
Einstein equations couple the UV and the IR, with small-scale fluctuations providing sources 
for the formation and evolution of large-scale perturbations. In this paper, we studied this UV- 
IR coupling of cosmological fluctuations. By integrating out short-wavelength perturbations we 
derived an effective theory for the long-wavelength universe. We observed that on sufficiently 
large scales the universe can be described by quasi-linear perturbations evolving in the presence 
of an effective fluid whose properties are determined by the small-scale structures. The fluid is 
somewhat unconventional in the following sense: higher moments in the Boltzmann hierarchy are 
small not because of sizable interactions like in a conventional fluid, but because they haven't 
had sufficient time to develop during one Hubble time In the absence of gravity 1~L~ l would 

go to infinity and a macroscopic effective fluid description would not be applicable at non-zero 
momentum. It is in this sense that we refer to the fluid as a 'gravitational fluid', highlighting 
the fundamental importance of gravity. On superhorizon scales, the only effect of small-scale 
structures is to renormalize the background density and pressure by terms of order the velocity 
dispersion. Moreover, we proved that the virial theorem naturally filters out the contributions 
from very small scales. On subhorizon scales, dissipative effects like viscosity are induced by the 
small-scale non-linearities. The imperfect fluid is described by a few parameters like the equation 
of state, the sound speed and a viscosity parameter. We proposed that this effective description 

29 It is straightforward to show that, even hypothetically, one cannot prepare a dense gas of black holes. For a 
large number N of black holes, the number density has to be smaller than 1/N 2 in Schwarzschild-radius units. For 
denser configurations, the whole gas collapses to form a big black-hole. 
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of the long-wavelength universe be used to formulate a well-defined alternative to conventional 
perturbation theory. At scales larger than some smoothing scale A -1 the theory is defined as an 
expansion in small perturbation variables - 5(k <C A) and 8(k <C A) - which evolve in a fluid 
whose physical parameters are determined by a numerical matching calculation. It remains to be 
quantified if a simple version of this procedure (e.g. a one- loop calculation in the effective theory) 
will be sufficient to reach the required level of precision in a practical application like baryon 
acoustic oscillations [60]. 
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A The Effective Fluid a la Newton 



In this appendix we give a simple derivation of the effective stress-energy tensor by smoothing 
the Newtonian conservation equations. 

A.l Preliminaries 

Consider a collisionless gas of non-relativistic particles. We define the single-particle phase space 
density f(x,p,t), such that fd 3 xd 3 p is the number of particles in an infinitesimal phase space 
volume element. Here, t, x are the time and spatial coordinates in the Newtonian frame defined 
in §2, and p = mx is the conjugate momentum. Phase-space conservation is described by the 
collisionless Boltzmann equation (or Vlasov equation) [23], 

^_V + I p .V/-mV».^-0, (A.l) 
Dt at m op 

where we have used the particle equation of motion 

* = -mV$. (A.2) 

Here, $ is the Newtonian potential defined as in §2. Taking moments of the distribution function 
/, we define the particle density, p m , the peculiar velocity flow, v, and the stress tensor, k^, 

ml f(x,p,r)) = f>m(x,t), (A. 3) 

f Pi 

m / —f(x,P,t) = p m Vi(x,t), (A. 4) 
Jp m 

/ -^f f(x,p,t) = p m ViVj(x,t) + Kij(x,t) . (A. 5) 
Jp m 



m 



Taking moments of the Boltzmann equation (A.l) gives equations relating the fields in Eqns. (A. 3)- 
(A.5). The zeroth moment gives the continuity equation (describing mass conservation) 

d t Pm + V i (p m v i ) = 0. (A.6) 

The first moment gives the Euler equation (describing momentum conservation) 

p m [dtv 1 + v>*] = -p m Vi$ - V jKij , (A.7) 

while the second moment gives an evolution equation for the stress tensor Kij. The hierarchy of 
moments of the Boltzmann equation is often truncated by either assuming that the stress tensor 
is negligible on large enough scales, ~ 0, or by postulating a phenomenological ansatz relating 
to the other fluid quantities, e.g. 

r — 2 I 

Kij = -pSij + rj ViVj + VjUi - -5ijV ■ v + C<% V • v , (A. 8) 

where p denotes the pressure and i] and ( are viscosity coefficients (see Appendix C). In this 
appendix we show that even starting with vanishing k^, integrating out short- wavelength fluc- 
tuations induces a non-zero effective stress tensor of the form of Eqn. (A. 8) at long-wavelengths. 
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A. 2 Integrating out Small Scales 

We are interested in the theory at scales k much larger than the scale of non-linearities q*. 
To define this effective long- wavelength theory we integrate out short- wavelength modes below a 
scale A <C starting from Eqns. (A. 6) and (A. 7). In real space, integrating out short-wavelength 
modes corresponds to a convolution of all fields with a window function W A , 

Pe = [p] A (x) = [ W A (\x - x'\) Pm (x') , (A.9) 

Jx' 

fa = [<t>] A (x) = I W A (\x-x'\)<t>(x'), (A.10) 

Jx' 

p e vi = [pmvMx) = [ W A (\x - x'\)p(x'y{x') . (A.ll) 

Jx' 

Essentially, this amounts to averaging the fields over domains of size A -1 . The smoothed velocity 
field v\ was defined such that piv\ represents the center of mass of the domain. In Fourier space, 
Eqns. (A.9)-(A.ll) become 

pi = W A {k) ■ Pk , (A.12) 
4 = W A (k) • <t> k , (A.13) 
(p f v\) k = W A {k)-(pv i ) k . (A.14) 

The window function is normalized such that 

W A (k = 0) = J d 3 xW A (x) = 1. (A.15) 

For concreteness we sometimes specialize to a Gaussian window function 



W A (k)=exp[-~] . (A.16) 

We split all fields into long and short modes 

Pm = pe + Ps, 4> = (pi + 4> s , v l = v\ + v\, (A.17) 

where in Fourier space the short-wavelength fluctuations are defined by 

p% = F A (k)-p k , (A.18) 
<t>% = F A {k)-<t> k , (A.19) 

(PmvDk = F A (k) ■ {pmV^H , (A.20) 

with F A (k) = 1 — W A (k). For notational simplicity we will often drop the index V denoting 
short-wavelength quantities except in cases where this could cause ambiguity. 

A. 3 Effective Stress from the Euler Equation 

To analyze the effects of small-scale non-linearities on the long-wavelength universe, we consider 
the Euler equation describing momentum conservation in a pressureless fluid 

p m [v* + uJ'VjV] + p m Vi$ = . (A.21) 
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We note that the Euler equation is linear in the Newtonian potential <J> and quadratic in fluid 
velocities v l . Notice that we never expand in density perturbations Sp m . The density p m and the 
Newtonian potential <3? are related by the Poisson equation 

V 2 $ = 4vrGp m . (A.22) 

Eqns. (A. 21) and (A.22) are valid on subhorizon scales and receive corrections for modes with 
wavelengths comparable to the horizon. This is made explicit in the more formal relativistic 
treatment of the same results in §3. 

We now show that Eqn. (A. 21) can be written in the form 



Pi 

where 



4 + vjVjvi + Pe V i <l>e = -V j [T j i \ A , (A.23) 



t*j = Pmvtv* - — • (A.24) 

As before [■■■]a denotes a spatial average over domains of size A" 1 and the labels £ and s denote 
long- wavelength smoothed fields and small- wavelength fluctuations, respectively. We will inter- 
pret Tij as the effective stress tensor induced by products of small-scale fluctuations. It sources 
the evolution of the long-wavelength perturbations. 

The proof of Eqn. (A.23) is instructive and shows that the result is valid up to higher-derivative 
corrections in the long-wavelength fields. 

Proof: 

Consider first the smoothing of the term p m v % + PmV^VjVi, 

* = f W A {\x-x'\)[p m d t v\x') + Pm v^x')d f v\x')\ . (A.25) 

J x' 

Using the continuity equation d t p m = —dj{p m v^) this may be written as 

★ = d t {p t v\)+ f W A v i d r (p m v j ) + f W A p m v^d y v i . (A.26) 

Jx' Jx' 

Integrating the second term by parts and using dj>W A = —djW A , we find 

* - dt{p t v\) + dj f W KPm v l vi . (A.27) 

J x' 

We split the fluid velocities v l into long- wavelength modes v\ and short-wavelength fluctuations v l SJ 

v i (x')=v}(x')+v i s (x'), (A.28) 
and expand v\{x') in a Taylor series around x, 

wj(aj') - v\(x) + djv\(x) ■ (x - x') j + ^d k djv\{x) ■ (x - x'y(x - x') k + ■■■ . (A.29) 

The derivative terms in Eqn. (A.29) will result in higher-derivative corrections to Eqn. (A.23). In partic- 
ular, cross-terms between v\ and v\ will vanish up to higher-derivative corrections. For concreteness, we 
will use a Gaussian window function, such that 

dj,W A = -djW A =A 2 (x- xy\V A , (A.30) 
d v d y W K = didjW A = -A 2 6ijW A + A 4 (x - x')\x - x') j W A , (A.31) 
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but our results will only slightly depend on this example. Using Eqns. (A. 28), (A. 29), (A. 30) and (A. 31) 
we may show from the definition of v\ that 

/ W APm v l s = -d k v\( W A (x - x') k p m - \d k 8iv\ [ W A (x-x') k (x-x') l Pm , (A.32) 

Jx' Jx' 1 Jx' 

= -^^- 2 pe JS- (A - 33) 
Using the above derivative expansions we find 

/ W A Pmt>V=/ W A p m v}vi + 2 f W APm v< 1 , l vl ) + f W APm v\vl, (A.34) 

J x' J x' Jx' Jx' 



where 



f W A p m v\vi = p t v\4 + C\ j and 2 f W A p m v<fvj> = C l 2 3 . (A.35) 

J x' J x' 

Here, we defined the higher-derivative terms 

CY = Pl d ^-2 9kP ^^ + p l vf^, (A.36) 
- 2 d -^v^ - PiV f-^f ■ (A.37) 

/ W A W V = Pt v\4 + [ W A Pm v\vl + C ij , (A.38) 

J x' Jx' 



Hence, we obtain 



where 



Vvl ■ Vvl 



Notice that many of the higher-derivative corrections in Eqns. (A.36) and (A.37) have canceled. Eqn. (A. 27) 
becomes 

* = d t { Pl v\) + d 3 { Pl v\vl) + djlpmviviU + djC ij . (A.40) 
Expanding the time derivative using the smoothed continuity equation, 

/ W A [d tPm + dvifav*)] = d tPe + di( Pe vi) = 0, (A.41) 

Jx' 

we find 

★ = Pt v\ + P tv\dp\ + dj\p m vivi\ A + djC ij . (A.42) 
Next, we consider smoothing of the term Pm di$>, 

o = f W A (\x-x'\) Pm {x')d v <S>(x'). (A.43) 

J x' 

As before, we define a split into smoothed long-wavelength modes and short-wavelength fluctuations, 

Pm (x') = Pe {x') + Ps {x') , (A.44) 

<f(x') = Mx')+<Ps(x'), (A.45) 

where 

and jw K 4>.*-\**. (A.46) 
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Smoothing of the Poisson equation implies 

V 2 </><? = 4nGp e and V 2 </> s = 4ttG p s . (A.47) 

As before, Eqn. (A. 46) enforces the vanishing of cross-terms in Eqn. (A. 43) up to higher-derivative cor- 
rections, so that we get 

o = Pi di<l>t+ ! W A p s di,<f> s +V\ (A.48) 



where 

vi s djA^kPe (A 49) 

stands for all higher-derivative corrections. As before, there have been cancellations of some higher- 
derivative terms. Using the Poisson equation, V 2 </> s = AttG p s , we get 

o = p t diil>t + J W ^£^ d ^s + Vl ■ ( A - 5 °) 
Integrating by parts and using dj>W\ = —djW\, this becomes 

• - «** + »> L - L * « ■ 

After using the fact 

dj>d s di>di'6 s 1„ f T „ di'd s d. 



f d r ^ a dj>di><j>s _1„ / 



this reduces to 



A 

72, 



Using the Poisson equation, V 2 ^ = 4irG pt, the higher-derivative term (A. 49) may be written as 

Z>* = V,£> ,J , where Z>' J = — ^ ' fc ' - ^ . A.54 

8irG ■ A 2 

Adding the results of Eqns. (A.42) (★) and (A.53) (o) gives 

Pi [*i + v foj v i\ + PtVifa = -Vj K] A - V J -[ T j ;] sa , (A.55) 

where 

rV = - 87 rG ' (A - 56) 

and 



[^■^^^^^ + ^•'11:^" ^ - (A.57) 

Up to the higher-derivative corrections in Eqn. (A. 57) wc have hence arrived at the desired answer, 
Eqn. (A.23). QED ■ 
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A. 4 Renormalized Pressure and Density 

From the stress tensor (A. 56) we may read of the effective pressure induced at order v 2 , 

3 PcS = [r\] A = [pvMU - !^|^ . (A.58) 

Similarly, the background density receives corrections at order v 2 . In §4 we computed this effective 
density p e g directly from the second-order Einstein equations: 

PeS = h(v s )p m }\ + \ [Pm<p S ]A , Pm = j(v s )e^ 3<i>s p , j(v s ) = 1 , (A.59) 

2 V 1 - V s 

or 

(A.60) 

A 

Alternatively, the corrections to the energy density could be found by considering the first post- 
Newtonian corrections to the continuity equation. Eqn. (A.60) describes kinetic and gravitational 
corrections to the total energy density taking into account the relativistic boost to relate the 
physical number density in the inertial frame where the fluid velocity is v to the number density 
in the fluid rest frame. The final answer, Eqn. (A.60), is so intuitive that it could have been 
guessed. In §2 we showed that this form of p c g is consistent with the conservation of the effective 
stress-energy tensor, d^T^ = 0. 



PcS 



i 1 2 1 
1 + 2^ + 2 ( 
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B The Effective Fluid a la Einstein 



Our discussion in most of this paper has been formulated in a fixed gauge. The Poisson gauge 
was convenient because of its controlled small-scale properties in the Newtonian limit. In §B.l we 
collect standard results from second-order perturbation theory in Poisson gauge. Second-order 
gauge transformations can be subtle since the combination of two short-wavelength transforma- 
tions can induce a long-wavelength change in the metric. We cite the explicit transformations in 
§B.2. Finally, in §B.3 we show that our interpretation of the long-wavelength effective fluid is 
limited to gauges with controlled small-scale behavior such as the Poisson gauge. We define the 
conditions that have to be satisfied in order to define the effective fluid in the long-wavelength 
limit. 

B.l Second-Order Equations in Poisson Gauge 
B.l.l Metric Perturbations 

The metric in Poisson gauge is given by 

ds 2 = a 2 ( V ) [-e^dr/ 2 + 2^dxMr/ + (e" 2 %- + XijJdxW] , (B.l) 
with xu = (i.e. the trace is absorbed in $) and 

= Xij,i = ® ■ ( B -2) 

The gauge condition (B.2) eliminates one scalar degree of freedom from goi and one scalar and 
one transverse vector degree of freedom from gij. Thus, uii is a transverse vector, while Xij is a 
transverse-traceless tensor. 

Perturbations are split into first- and second-order terms 

* = + 2* (2) > * = + ^ > <* = "i 2) . XiJ = X ( 5 ) - (B.3) 

We have ignored first-order vector and tensor perturbations as they aren't produced by inflation 
and in any case decay with the expansion of the universe. The metric determinant is \J—g = 
e *-3<i> a 4 to seconc [ order. 

Christoffel symbols. The second-order Christoffel symbols are 



1 00 

1 0i 

pi 
1 00 

r°. 


= n + y, 
= y^ + nui, 


1, . 


(B.4) 
(B.5) 
(B.6) 

(B.7) 


-pi 


= (H- tySij + ]ptij + - u)j,i) , 




(B.8) 


r i 
1 jk 


fc (5j - $j<5| + $% k - uJb 3k + 


^ [Xij),k + (Xik),j ~ (Xjk),i ■ 


(B.9) 
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Second-order Einstein tensor. The Einstein tensor up to second order has the following com- 
ponents 

„-2* 



„2<I> 



= 2- 



+ (H — 



„2<J>+2* 



(B.10) 
(B.ll) 



(-{M 2 + 2W) - 2$^ - 3$ 2 + 2H(2<S> + *) + 2$) 



+ e 2 *(* fe * fc + - $ fcfc ) 



5 ) 



„2<I> 



1 

2^2 



+ 



1 

2^2 



xj + ^xj-Cxj),** • (B.12) 



We split these tensor components into G^ v = G^ u + (G M i,) L + (G M j,) NL . The background 
components are 

— fl 2 G^o = 3% 2 , = , 

The linear parts are 



-a't&j =n 2 + 2n. 



(B.13) 



-y(G u f = V 2 $-3^($ + ^), 



(n 2 + 2?i)* + n(2$ + + Jv 2 ($ - *)1 5) + 



4 L 



(wV+wf i ) + 2W( W V + 



or. 



+ 



x*. + 2Hx5--(xi- 



Finally, the quadratic parts are 
-a 2 (G° ) NL = 

2$[$ + fttf] ii -$tf i 



+ y(G*o) NL 



(B.14) 
(B.15) 

(B.16) 

(B.17) 
(B.18) 



+a 2 (G i 



\NL 



-4(n 2 + 2U)^ 2 - 2$tf - 3<£ 2 - 4^(2$ + - 4tf$ 
+ fc tf fc -2$[$-tf] )fcfc ]$j 



(B.19) 



In the absence of anisotropic stress c/> = $> = * and using c/)^ 3> %c/> ~ c/> on small scales we 



find 



-a 2 (G° ) NL « -^ + 4# # , 
+a 2 (G^.) NL « MfiSi-tyrfj- 



(B.20) 
(B.21) 



Notice that in this gradient expansion Eqns. (B.20) and (B.21) hold independent of the back- 
ground equation of state. 
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B.1.2 Matter Perturbations 



We now cite results for the parameterization of perturbations in the dark matter energy and 
momentum. We first give results for dark matter in the fluid approximation and then present 
the energy-momentum tensor for point particles. 

Dark matter in the fluid approximation. We consider the stress-energy tensor of a fluid in the 
covariant form 

T% = ( P + p)u»u v + p 6£ + £"„ , (B.22) 

where g^u^u 11 = — 1 and 

= = . (B.23) 

Eqn. (B.23) implies that only the spatial components of Yf v are non-zero, i.e. S° = S l = 0. In 
this paper we will often restrict to a background without pressure and anisotropic stress, i.e. p = 
E l ■ = 0, but our results are easily extended to more general backgrounds. At second order, 
fluid velocities and gradients of the gravitational potential induce both pressure and anisotropic 
stress. A characterization of the second-order corrections to the effective fluid is one of the main 
objectives of this paper. 

To second order in metric perturbations the four-velocity is 

u ° = a - 1 e-* 7 (v) , u l = a^eV , (B.24) 
uq = —ae^j(v) , Ui = a(uJi + e _(J V) , (B.25) 



where we defined the Lorentz factor 



7(«) = ^l + -v 2 . (B.26) 



Eqn. (B.22) becomes 



A = -J 2 p = -p(l + v 2 ), (B.27) 

r = -e*+V, (B.28) 

T° = p(ui + e-*-*«i) , (B.29) 

TV = pv^j. (B.30) 

Notice that the conservation equation, V^T M1/ = 0, may be written as 

u v V^ v = u u V^{pu»u u ) = -V^pu") = -d^yf^jpui*) = . (B.31) 

This motivates introducing 

Pm - = l(v)e-™p « p (l + £ - 3$) , (B.32) 

such that Eqn. (B.31) implies 

— \ d 3 xa 3 p m = 0. (B.33) 
drj J 
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The physical role of p m becomes clear when thinking about dark matter in terms of particles 
and their conserved number density, n oc p m . A relativistic boost relates the physical number 
density in the intertial frame where the fluid velocity is v to the number density in the fluid rest 
frame, n rcst oc p. The factor 7(v)e -3 * in Eqn. (B.32) clearly relates the volume element in the 
rest frame to the physical volume in the moving frame. 

The conserved energy of a localized system is (c/. Exercise 20.5 in Misner, Thorne and 
Wheeler [104]) 



p(l-v 2 



-1/2 



1 



1 



Lorcntz 



- - I dx 3 a 3 



J- n ~~ 



+ 2^ + 2 P ® 

: 

(G u ) 



KE PE 
\NL- 



(9xx9yy9zz) 1/2 dxdydz 

v ' 

proper volume 



8vrG 



/■ 

dx 3 a 3 Pm (l + -v 2 + -$) . 

More generally, the conserved four-momentum of a localized system is 

(^ ) NL - 



= - J dx 3 a 3 



8ttG 



d 3 x a 3 r M n 



(B.34) 

(B.35) 
(B.36) 

(B.37) 



For an early discussion of conserved energy and momentum in General Relativity see [105, 106]. 

Dark matter as point particles. When relating our approach to iV-body simulations a particle 
description of dark matter is essential. A collection of particles 'a' is described by the following 
stress-energy tensor [18], 

utul 5 D (x - x a (rj)) 



rn 



a n 



To second order, its components are 



1 o 



-a" 3 e 3 * 



E 



'-9 



1 + y J S D (x ~ x & (r})) , 



r = a- 3 e 3 *e*+*J2 m -< S »( x - x ^)), 



3„3* 



(B.38) 

(B.39) 
(B.40) 
(B.41) 



The spatial integral of the zero-zero component of the effective stress-energy tensor gives the 
total energy of the particles 

, (0°») NL ' 



Es -f 



l 3 xa 3 T° 



d 3 xa 3 



E 



u 8ttG 



Particles evolve according to the geodesic equation, which at second order is 



d OS pt , _ , ■ • , dtXj 

— ^ + {u-2<s>-m)— 



drj 1 



drj 



(B.42) 
(B.43) 

(B.44) 
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B.1.3 Einstein Equations 

Metric and matter perturbations are related by the Einstein equations 
In the following we present the first-order and second-order equations. 



(B.45) 



First-order solutions. In linear perturbation theory, the Einstein tensor is linearized in the 
metric perturbations <p = and ip = vt^ 1 ) 



V 2 4>-3H(<f> + nij) 
2. 



-4vrGa 2 (T° - f° ) 



4vrGa 2 T Q , 



didj 



1 



8vrGa 2 didj 



1 Hrpk 



where we have subtracted the background 

8vrGa 2 ^ n 



i 



U 2 + 2U 



SvrGa 2 -, 
3 *■ 



(B.46) 
(B.47) 

(B.48) 

(B.49) 

(B.50) 



In addition, the energy-momentum tensor may be linearized in the fluid three- velocities v % and the 
density contrast 6. The Einstein equations then simplify considerably: e.g. for a background with 
constant equation of state w = p/p, and in the absence of anisotropic stress (Sjj = => <f> = ip) 
the Newtonian potential obeys the following evolution equation 



4> + 3W(1 + w)<f> - wV 2 (f> = , 



(B.51) 



where we have used Eqns. (B.46) and (B.48) and 5p = c 2 Sp with c 2 = w. In Fourier space 
Eqn. (B.51) becomes 

6(1 + w) 1 : . l2 



1 + 3w rj 



4>b + wk 0k = , 



which has the following exact solution 



MV) = V a [Ci(k)J a (y) + C 2 (k)Y a (y)} , y = ^nkn , a = 



1 fb + 3w 



2 V 1 + 3w 



(B.52) 



(B.53) 



where J a and Y a are Bessel functions of order a. During the matter-dominated era, w = 0, this 
becomes 

C 2 (k) 



whereas during the radiation-dominated era, w = \ , we find 



y 



c, 



(*)( 



smy 



cos y 



)+C 2 (fe)( 



cos y 
1/ 



+ siny 



(B.54) 



(B.55) 



In both cases the decaying mode may be dropped by setting C^A:) = 0. For a radiation- 
dominated background, the Newtonian potential is time-independent on superhorizon scales, 
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lim/j^o <t>\i(v) = Ci(k), but decays on subhorizon scales. In contrast, during matter-domination 
the growing mode linear gravitational potential is time-independent on all scales, with a spatial 
profile 4>(x) given by the Poisson equation 

V 2 4> - m 2 <j) = AnGa 2 p5 , (B.56) 



or 

-k 2 (f>k = \u 2 5 k + m 2 (t> k . (B.57) 
The linear velocity is determined by the gradient of the potential, see Eqn. (B.47), 

v i = ~^i- (B.58) 

Using the scalar velocity potential, Vi = diV, Eqns. (B.56) and (B.58) may be combined into 

-k 2 fa = \n 2 d k , (B.59) 

where 

d k = S k - mv k (B.60) 

is the gauge-invariant linear density perturbation in the fluid rest frame. Eqn. (B.59) is valid on 
all scales, while the Newtonian approximation, —k 2 (f) k ps ^T~L 2 5 k , is only valid for scales that are 
much smaller than the horizon. 

Second-order theory. The non-linear Einstein equations can be written in a form that is very 
similar to the linear Eqns. (B.46)-(B.49) if we replace the stress-energy tensor Tt by the stress- 
energy pseudo-tensor rt: 

n K = Tt - [ -^L_ . (B.61) 

The Einstein equations for $ = + 2$ (2) and = + ^ (2) then become, 

(G^) L = 87rG(^-ft), (B.62) 

or explicitly, 

V 2 $-3^(6 + W) = -4vrGa 2 (r° - f° ) , (B.63) 
[6 + W&] . = AnGa 2 r\ , (B.64) 

$ + ft(2$ + *) + (ft 2 + 2?i)tf- ^V 2 ($-*) = ^^-{t\-t\), (B.65) 



didj 



($-n»i-3^($-*) 



= SnGcPdidj 



1 



■i k 



(B.66) 



This defines the dynamics of long-wavelength scalar fluctuations sourced by products of short- 
wavelength fluctuations. All gravitational non-linearities have been moved to the r.h.s. of the 
Einstein equations. Eqns. (B.63) and (B.64) can be combined into the Poisson equation, 

V 2 $ = 4^Ga 2 [mV^diT^ - (r° - f° )] . (B.67) 
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At long wavelengths, rt is conserved by virtue of the linearized Bianchi identity 



(B.68) 



Using (G^) L = 87rG(r^ — rt) and computing the linear covariant derivative action on the 
background this is 



-p 



3$ <L o + £„ 



(B.69) 



Notice that if one wished to obtain this equation from the Newtonian stress tensor we found in 
§2, it would be necessary to keep subleading terms in x l in the expression for t^ v . 



B.1.4 Conservation Equations 

Finally, we give the velocity expansion of the continuity equations, V^Tt = 0. At order v 2 the 
mass conservation equation and the Euler equation are 

\^-{a 3 Pm) + -Viipmv') = , (B.70) 

and 

— (avi) + v j V jVi = -Vi* , (B.71) 

where p m was defined in Eqn. (B.32). A common way to organize the post-Newtonian expansion 
around Minkowski space is as an expansion of the field equations in powers of c _1 (see e.g. Chan- 
drasekhar [57]). For cosmology this post-Newtonian expansion at order 1PN (c~ 2 ) was performed 
in [107] . Our velocity expansion at order v 2 effectively keeps 1PN order in the continuity equation, 
but Newtonian OPN in the Euler equation. 

B.2 Second-Order Gauge Transformations 

To relate the results obtained in Poisson gauge to other gauges we need to consider second- 
order gauge transformations. As a reference for the reader we here cite those transformations for 
metric perturbations. More details may be found in the fantastically useful review by Malik and 
Wands [52]. 

The most general set of metric perturbations can be written as 

8 goo = -2a 2 ifj , 5g 0i = a 2 Bi , 5gij = 2a 2 Cij , (B.72) 

where 

Bi = B,i + Ui , (B.73) 
Cij = -4>Sij + E tij + F (lJ) + ^Xij ■ (B.74) 

Here, ip, B, (ft, and E are scalar perturbations, LOi and Fi are vector perturbations, and Xij 1S a 
tensor perturbation. In this parameterization the Poisson gauge is defined by B = E = Fi = 0. 
Consider the second-order coordinate transformation 

^ = ^ + % + ^i)^ri)+%)- (B.75) 
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At every order the four-vectors (r = 1, 2), can be decomposed into scalar and vector parts 
#)=«(,■), ^ r) = d*/3 (r) + ( f (r) , with did{ r) =0. (B.76) 

Second-order gauge transformations of tensors are written in a compact form by the use of Lie 
derivatives. Any tensor (e.g. the metric g) transforms as 

~g = g + C kl) g+ l -{Cl m +C k2) )g, (B.77) 

where we have suppressed indices. 30 Using 



this implies 



9 = 9 + Sg w + -Sg {2 ) , (B.78) 



<%) = $9(1) (B-79) 
*5(2) = 8g ( 2 ) +£t m g + £l (1) g + 2£z w 6g w . (B.80) 



B.2.1 First-Order Transformations 

We find the following first-order transformations: 
• scalar perturbations 



• vector perturbations 



tensor perturbations 



fa = 0i -Uai, (B.81) 

^ = 1 p l+ '}{a l + a l , (B.82) 

Bi = B 1 -a 1 +p 1 , (B.83) 

Si = Si+/3i, (B.84) 



= wj-dj, (B.85) 
F{ = Fl + d\, (B.86) 



Xiij = Xiij- (B.87) 



Furthermore, the scalar shear potential, o\ = E\ — B±, and the momentum scalar, q\ = v\ + i?i 
transform as 

o-i = <7i + ai , (B.88) 
?i = gi -ai. (B.89) 



'Explicitly, the Lie derivative of the metric tensor is 



2 , 



From this follows straightforwardly an expression for C^g^ v = C^C^g^v) 
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B.2.2 Second-Order Transformations 



Second-order gauge transformations are considerably more involved since they receive 
tions from products of first-order transformations. 
We define 



= OL\ 



on + bUai + (H + 2H)ai + m^i + 2ipi 



+2di(di + 2^i) + £i fc (di + Uai + 2^i), fc + ^ ai, fc - 2B ifc - £ lfc 



0i = 2 [(2WB M + B u )ai + B u ^\ - 2^ia M + B lfc ^ + B u di + 2CWJ 8 
+4?tei(£ii - ai,i) + - 3ai,j) + - di,j) 

+H + 2 £ifc,j) + £1 (£ii,Jfc - «i,ifc) - ai,fc£i,i , 



and 



= 2 



(H 2 + 2^)a? + W ( aidl + aik ^) 



Sij 



-2a hi a ld + 4ai(C , Hi + 2HC Uj ) + 4 J B 1(i o; J) 

+ 2 «i£i(ij) + 2 £i(i,j)fc£i + 2 £i(i,fe£ij) + 2 £i(i«ij) 
The second-order metric transformations then can be written as follows: 
• scalar perturbations 



92 = 

V>2 = 

B 2 = 

Eo = 



4>2 — I-L012 



4 " 4 
ip2 + 7-Lot2 + d 2 + E , 
£ 2 - a 2 + /3 2 + V~ 2 6 fc fc , 

e 2 + h + ^v- 2 v- 2 r j ' . . 



_ v -2 T fc 



k > 



vector perturbations 



w 2i 



w 2i 



2i 



Qi + v- 2 e fc 



2^ 



= Fa + <fe + V^T^ - V~ Z V~ Z T 



• tensor perturbations 

X2v = X2^ + Ty + i(v- 2 T%-T\)^ + ^V- 2 V- 2 T w 



+^V- 2 T^ 



V" 
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In the case of a pure time-shift, & = 0, Eqns. (B.90), (B.91) and (B.92) become 



&i + 5Ha! + (% + 2U)a x + 4^^i + 2^i + 2Ai(di + 2^i) , 



e, = 2 



{2'HBii + B u )ai - 2^iai,i + B Xi a x 



4Ha\ai i — ?>c\\a.\ i — ct\o.\ 



and 



Tij — 2 



(H 2 + 2W)a? + Uaxax <5y - 2ai,iaij + 4ai(C Hj + 2^C W ) + 4B 1(i a J) . (B.102) 



(B.100) 



(B.101) 



B.3 Covariance of the Effective Theory 

Next, we discuss the gauge conditions that have to be satisfied in order to define the effective 
fluid in the long- wavelength limit. 



B.3.1 Long- Wavelength Limit of the Metric 

Let us split the metric into a background g^ u (k — > 0), long- wavelength perturbations H^ik <C A) 
and short-wavelength perturbations h^ v {k 3> A), i.e. 



(B.103) 



The long-wavelength limit, k — > 0, extracts the perturbations and projects out the 
perturbations 



g^ = lim g^ u = g^ u + H^ u . 

fc<A 



At linear order in Hn U and quadratic order in hn U we find 



-9 



1 



1 



1 



g,- = ^ -h^ -bT + h> ta h v a , 
which in the long-wavelength limit becomes 



lim • 



-9 



-q I 1 + -H + lim 

1 2 fc«A 



\h 2 ~ -fihi 



lim g^ 



(B.104) 

(B.105) 
(B.106) 

(B.107) 
(B.108) 



Notice that the long- wavelength limit of the covariant metric component g^, Eqn. (B.104), is not 
the inverse of the long-wavelength limit of the contravariant metric component g^ u , Eqn. (B.108), 
unless the long-wavelength limit of the product of short-wavelength perturbations is sufficiently 
small. In other words, to define a covariant effective theory in the long-wavelength limit we require 
the following condition to hold 



lim \H^, \ » lim h^h^R g a e , 
fc«A' ( 2 )' fe«A W (l) yap ' 



(B.109) 



i.e. it needs to be the case that the k — > limit of the second-order piece, H( 2 ) ls bigger than the 
k — > limit of the square of the first-order piece, h 2 ^. In this case the long- wavelength metric is 
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well-defined and can be used to construct covariant actions such as the action of a scalar field tp 
living in this effective spacetime g^ u , 



-ggrdpipdw. (B.110) 
We will now discuss in which gauges the condition (B.109) is satisfied. 

B.3.2 Poisson Gauge 

Let us first demonstrate that Poisson gauge is a 'good gauge' in the sense defined above. In 
Poisson gauge the scalar perturbations of the metric are 



<$goo = -2a 



Sgij = -2a 2 



4 1 ' + \^ 



(B.lll) 



Sij , (B.112) 



5g 0l = 0, (B.113) 
where during pure matter-dominance we find (see also Refs. [2, 11]) 

*W = $W = < j ) = const j (B.114) 

and 



,(2) _ 2 1 o , Cl *(2) _ 2 1 

Here, we have defined 



*p =*5^ + C*> ^ = 5 + ^. (B.115) 



S = * V~ 2 [2(c^) 2 + 5(V 2 0) 2 + 7drfdiV 2 <p] , (B.116) 

and the constants Cy, and are set by initial conditions [2, 11]. For the following estimates we 
will use 

(B,17, 

It is then easy to check that the condition (B.109) implies 

\di<p\ 2 > H 2 (f> 2 . (B.118) 

Eqn. (B.118) is satisfied if we couple two subhorizon perturbations with k > %. This shows that 
integrating out subhorizon modes in Poisson gauge gives a well-defined long-wavelength metric. 

B.3.3 Arbitrary Gauges 

Of course, we don't expect the condition (B.109) to hold in all gauges. We therefore now exhibit 
second-order gauge transformations of our Poisson gauge results to define a set of 'good gauges' 
in which a consistent effective fluid description is possible. 

Second-order gauge transformations. Recall the gauge transformations of the metric tensor 

<%) = Sg w +C ({1) g, (B.119) 
5g { 2) = Sg (2) +C k2) g + Cl w g + 2C kl) 5g {1) . (B.120) 
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The most general second-order gauge transformations were presented in §B.2. Here, we will 
consider the special case of a pure time-shift, C = 0, and focus on the £f {1) <? term in Eqn. (B.120) 
(inspection of the results of §B.2 shows that this captures the leading effect in an expansion in 
scale-scale gradients). For the FRW background 

goo = -a 2 , g i = , = a 2 <% , (B.121) 

we may compute £ 5(r) <? and £| 0, using C^g^ = g^ u ,aC a + £^<W + £^<W We find 

£$900 = -2a 2 [a + Ha] , 



£$90i = -a a,i , 
C^cjij = +2a 2 Ha , 



and 



H ( i)50o 

^£ (1) 0Oi 



= -2a 2 



= —a 



= -2a z 



a[a + 5Ha + {% + -H 2 )a] + 2d 2 
-a[(?i + % 2 )a + Ha] + ^(a ;) 2 



Sij + • • 



(B.122) 
(B.123) 
(B.124) 



(B.125) 
(B.126) 
(B.127) 



We will be interested in first-order time shifts on very small (subhorizon) scales. In that case 
the dominant terms are the ones with the maximal number of spatial gradients. In this gradient 
expansion we find 



and 



£(9oo = a 2 0(7La) , C^g i = -a 2 a ti , C^gij = a 2 0(7La) 



C( (1) goo = a 2 (D(H 2 a 2 ), C^fa = a 2 0(Haa,i) 
1 



r 2 - 



-2a z 



-Xatf + 0{U 2 a 2 " 



Sh + 



(B.128) 

(B.129) 
(B.130) 



Second-order time-shift and 'good gauges'. Let us now go from Poisson gauge to an arbitrary 
gauge using the second-order gauge transformations discussed above. 
We find, e.g. 



6g$ C -a 2 a. 



and 



^C-2a 2 [4 2 ) + i(a^ 



Consider then 



and 



The condition 



\H {2) \ (Z5~a {2) 



"J 



1(1)5,(1) gap c _ a -2l(l)l (1) _ Q 2 a . 



H (2) \^>h (1) h (1) a afi 



ia j/3 



(B.131) 
(B.132) 

(B.133) 
(B.134) 
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now implies 

$?»(a,,) 2 . (B.135) 

Only gauges that are connected to Poisson gauge in this way have a well-defined long-wavelength 
limit for the metric after integrating out short-wavelength modes. This means that we only reach 
a 'good gauge' from Poisson gauge if |9ja| 2 does not give the dominant effect in the new metric 
perturbation, i.e. 7i 2 \dia\ 2 <C |9j0| 2 . A popular gauge that violates the condition (B.135) is the 
comoving gauge, T* = 0. In that case, %a ~ <fi, and 

\H\f\~h { Sh { $9 aP - (B.136) 

In comoving gauge we therefore cannot define a unique long-wavelength metric after integrating 
our short-wavelength fluctuations. This is related to the breakdown of that gauge on small-scales. 
Analogous considerations apply to the second-order metric in uniform density gauge. 

B.3.4 Background Field Method 

An elegant alternative to ensure covariance of the theory in the long- wavelength limit is the back- 
ground field method (BFM). The background field method is a technique for quantizing gauge 
theories without losing explicit gauge invariance of the effective action (see e.g. Ref. [108]). The 
fields of the classical Lagrangian C are decomposed into background fields (p (analogous to g) and 
quantum fields <p (analogous to h). While the background fields are treated as external sources, 
only the quantum fields are variables of integration in the functional integral. A gauge-fixing 
term is added which breaks only the invariance with respect to quantum gauge transformations, 
but retains the invariance of the functional integral with respect to background field gauge trans- 
formations. From the functional integral an effective action S'cfrf^] for the background fields is 
derived, which by construction is invariant under gauge transformations of the background fields 
and thus gauge- invariant. In this formalism the effective action is manifestly covariant after 
integrating out the short-wavelength modes. 

The BFM has been applied to perturbative General Relativity by Feynman, DeWitt, 't Hooft 
and Veltman, e.g. Ref. [109]. As before, the metric is split into a smooth background, g^v, and 
short- wavelength perturbations, h^ u . Ultimately, we want to derive the effective action for g^ u 
by integrating out quadratic fluctuations in h^ u . The key element of the BFM is a gauge-fixing 
term that is covariant with respect to the long- wavelength metric g^ u , 

£ g augc = f M f 11 , f^t = D u ( V - , (B.137) 

where is the covariant derivative with respect to the background metric g^. This ensures that 
the effective action for g^ v is covariant. It may be shown that at leading order in our gradient 
expansion the BFM with the gauge- fixing Lagrangian (B.137) is equivalent to working in Poisson 
gauge, (f M ) P ps 0. 
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C Dissipative Fluid Dynamics 



The theory of imperfect fluids is reviewed comprehensively in the books by Weinberg [18] or 
Landau and Lifshitz [58]. For the benefit of the reader we here collect the main elements of that 
analysis. In §5 we showed how this connects to our theory of the gravitational fluid. 
For imperfect fluids the stress tensor may be written as 

Tij = pUiUj + {p- (0)-Yij + Sjj , (C.l) 

where is the viscous stress tensor describing dissipative dynamics. At leading order in an 
expansion in spatial derivatives an ansatz for E^- that is consistent with the second law of ther- 
modynamics is [18] 

= -Wij , (C2) 

where r] > and ( > are the coefficients of shear and bulk viscosity, respectively and 
(7^ = — \5ijVk,k- With Eqn. (C.2) the Euler equation becomes the non-relativistic Navier- 
Stokes equation. In general, for fluids that are locally at thermal equilibrium one also has terms 
proportional to the temperature gradient, weighted by the heat conduction. Though we expect 
to have dissipative terms also in tqi that can in principle correspond to temperature gradients, 
our effective fluid is clearly not in thermal equilibrium, so it is not clear what these terms would 
correspond to in our case. Fortunately, we do not need such an interpretation: as we saw in the 
main text, to describe dissipative effects in our case, we just need the higher-derivative expansion 
of the stress tensor — cf. Eqn. (5.2). 

In an adiabatic fluid the pressure is a unique function of the energy density, p = p(p), and 
the adiabatic sound speed is fully determined by the equation of state w = p/p, 

C n — - = W — r . (C3 

a p 3H(l + w) v ' 

In general, the pressure might not be a unique function of the energy density. The sound speed 
is then defined as the ratio of pressure and density perturbations in the frame comoving with the 
fluid 

The speed of sound may be extracted by correlating the pressure with the density contrast of a 
long- wavelength perturbation 5i, 

We define the scalar component a of the anisotropic stress tensor Sjj as 

(p + P)<J = V- 2 (didj - ^ tJ V 2 j . (C.6) 
We may relate a to the shear viscosity in Eqn. (C.2), 
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where we introduced the tilde viscous speed of sound, c v i s (as for c v j s , we stress that this is a 
dissipative term that does not induce propagating sound waves). We see that the parameter 
c v j s rnay be determined by correlating the anisotropic stress with the velocity divergence of a 
long- wavelength perturbation 6g, 

^ n Wi)- (c - 8) 

For c vi > this adds a dissipative term to Euler equation. 
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D Estimates in Perturbation Theory 



In this appendix we illustrate some of the ideas of this paper by presenting sample calculations 
in tree-level perturbation theory. In particular, we wish to show how the parameters of the long- 
wavelength effective theory can be computed if we know the theory that describes the dynamics 
of the short-wavelength universe. This is the standard procedure of matching the parameters 
between two effective theories that are valid at different scales. In this appendix, as an example, 
we will implement this procedure by assuming that perturbation theory is applicable to describe 
the small-scale dynamics in a dark matter universe. 

D.l Fluid Parameters at Tree Level 



Example: 

Perturbation theory. During the matter era, the method of perturbation theory is essentially to perform 
an expansion in the scale factor a(rj) [24] 31 



5(k, v ) = J2a n ( V )S n (k), (D.l) 

71=1 

CO 

0(k, V ) = -n( V )J2 an (v)0n(k), (D.2) 



n=l 



where S n and 9 n are time-independent mode-coupling integrals over n powers of the initial density field; 
e.g. at second order, 

S 2 (k)= f f S B (k-q 1 -q 2 )F 2 (q 1 ,q 2 )S 1 (q 1 )5 1 (q 2 ), (D.3) 

Jq 1 Jq 2 



F 2{qi , q2 )^ 1 + \^ + f) + 2 ^, (D.4) 



where 

5 , 1 fqi o 2 \ . 2 2 q l ■ q 2 

The angular integral of the kernel results in a pure number 

V2 ^\ ^r F *(<ii,<i 2 ) = //^ = ^. (d.5) 

We furthermore define 

ee afa^fc) , <J< 2 >(„) ee a\r,)5 2 (k) , (D.6) 



such that 



4 = 4 1) +4 2) - (D.7) 



The convolution in Eqn. (D.3) captures the UV-IR coupling between long-wavelength modes and short- 
wavelength fluctuations. Using these results we may evaluate the correlation functions in the definitions 

of p eff , c s and c vis , 

1 ,. / \ 2 1 (SePes) ~ 2 (<Weff) m c x 

Peff = 3 I"?o< t) ' Cs iw Cvis = l^P (D ' 8) 

where 



: In a general cosmology, the scale factor expansion is replaced by an expansion in the growth factor D + (a). 
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These formulas assume negligible bulk viscosity, which otherwise would contribute to the trace of the 
stress tensor. As we remarked in the main text, at tree level in perturbation theory it is impossible to 
distinguish between bulk and shear viscosity. In this appendix we only consider tree-level computations 
because our purpose is simply to explain how the various terms of the long-wavelength effective theory 
can be extracted from perturbation theory in principle. 

Renormalization of the background. We first consider the renormalization of the background (see §4), 
i.e. the fluid in the limit of long wavelengths. In perturbation theory we may use the linear relation 
between the fluid three- velocity and the Newtonian potential, Eqn. (B.58), to write 



Wa£g&- (DA0) 



and 



From this we find 



8nGa 2 ' 1 3 8irGa 2 ' 



(D.ll) 



M£ = -W A (p)-^ I a(q,p)S q 6 p _ q , a{q,p)= q ^_^ 2 , (D.12) 



and 

/2 



(Hp) = -W A (p) (2n) 3 S B (p) Vl± f a(q,p)P s (q) . (D.13) 

Jq 

The effective pressure of the background is 32 

PcS = 1 Jm([r]*) = £ jf dlngAfo), (D.14) 
where we defined the spectral density of velocity fluctuations (see Fig. 3) 

A 2 (g) = ^A 2 (g)^A 2 (g). (D.15) 



This shows that the renormalization of the background pressure is always positive and scales with the 
velocity dispersion of small-scale fluctuations. We remind the reader that the integral in Eqn. (D.14) is 
restricted to short-wavelength modes, q > A, or 

J^d\nqA 2 v {q) = J din q I%(q) A 2 v (q) , (D.16) 

where F A {q) = 1 - W A {q). 

We extract the scalar component of the anisotropic stress as follows 

°M = \ P -f([r%), (D.17) 

where 



if we define 



^in-]p = -W A (P) ■ \n 2 p f p(q,p) 8 q 6 p _ q , (D.18) 
P Jq 

P(q,P) = YY~t — ^2 • D. 19) 

p z q z (p — qy 



32 Eqn. (D.14) assumes negligible bulk viscosity which otherwise would also contribute to the trace of Tij. As 
we explained in §5 the presence of significant bulk viscosity requires going beyond tree- level perturbation theory 
and/or numerical simulations to measure the individual viscosity parameters separately. 
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Figure 3: Spectral density of velocity fluctuations, A^(g), computed in linear perturbation theory. Non- 
linear corrections are significant in the regime q > <7nl- The spectral density A^(q) appears in the kernel 
of all convolution integrals considered in this paper. 



Hence, we find 
Using 

this implies 



-W A (p) (2^) 3 5 D (p)-H 2 I P{q, P )Ps{q) 



lim iB(q,p) 

p-i-0 



cr off = lim cr G ff(p) oc 

p— >0 



ll 



(D.20) 
(D.21) 
(D.22) 



This confirms our arguments in §4 that the effective anisotropic stress vanishes on very large scales. 

Tidal pressure. We now consider the effects of long-wavelength perturbations on the properties of the 
effective fluid. To this end, we compute the correlation of the effective pressure with a long-wavelength 
mode 5 e k | k < A, 



(Sep, 



1 



'cflV 



I r_iA\ 
pi 



(Si [r 



-W A {p) 



12 



a{q,p) (S e k S s q S s p _ q ) 



(D.23) 



The correlation function between a long-wavelength mode, Sj, = #1 , and the product of two short- 
wavelength modes, 5* = Sq 1 ^ + S q 2 \ is 



(sls q s;_ q ) = <4 1) 4 2) 4 1 -,) + (4 1) 4 1) 4 2 - g ) 
= 2x<4 1 ^>4 1 ) q ), 



(D.24) 

■„ •>-,,- (D-25) 

where we have used the fact that correlation functions with an odd number of 6^ vanish if the initial 
fluctuations are Gaussian. Since S q ' is written in Eqn. (D.3) as a convolution of a long-wavelength mode 
and a short-wavelength mode, the three-point function in Eqn. (D.25) is related to a four-point function 
of first-order fluctuations, 



where 



(,5 (1) <5«5«<5 (1) ) 
\ ,fe ,9i °fl2 p -q> 



Sv(q -qi~ q 2 )F 2 ( qi ,q 2 ) (6^S^S^S^_ 
2 x 



(D.26) 

(D.27) 
(D.28) 
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Using 

= (2nfP s (k)6 D (k + gi ), 

= (2vr) 3 P 5 (|p-g|)^(p-g + 9 2 ), 

we perform the q 1 and q 2 integrals in Eqn. (D.26), 

^ 2) 4-«) = ^f p s{k) + p) 2 (2^) 3 F 2 (fc, fe + 9) Ps(\k + q\) . 



(D.29) 
(D.30) 

(D.31) 



We find 



, s , -W A (k)- { -^n 2 p^P S (\k + q\)a(q,-k)F 2 (k,k + q). (D.32) 

Using that lim^A Wa(&) ~ 1 an d that the integral is dominated by modes with q^> k, this simplifies to 

(S ePcS ) (2tt) 3 r 



{StSt) 



-p f P S (q)%F 2 {k iq ). 
J a Q 



or 



2 _ 1 (jl_geff) _ V2 
Cs ~ p (6 t S t ) ~ 3 



dlnqA 2 v {q) 



(D.33) 



(D.34) 



Anisotropic stress and viscosity. Next, we compute the correlation of the anisotropic stress with a 
long- wavelength mode Sf, = —9 i k /'H. \ k < A, 



a cff ) = ^<#[tV]£) - -W A (p) • / , (D.35) 

p P A Jq 



where /3(q,p) was defined in Eqn. (D.19) and the three-point function was computed above, 



(SU'6') = (SiS") 4 (27r) 3 F 2 (fe, q)Ps(q) . 



We find, 



(tie <r e g) 



= -10(2ir) 3 H 2 I Ps(q)F 2 (k,q)f3(q,-k). 

J n 



(S t 6 t ) 

As before, this simplifies since modes with g> k dominate the integral 



(St <7 e ff) 

(StS e ) 



{k-q) 2 -\k 2 q 2 



= 10(2tt)^ / P 5 (q)F 2 (k,q) 



or 



where 



^ = %# = 10^jfdln^( ff ) 



M 2 = |^F 2 (fc,q)((fe.g) 2 -i) 



1\ 16 
315 



(D.36) 
(D.37) 

(D.38) 
(D.39) 

(D.40) 



The ratio of the sound speed associated with tidal pressure and the sound speed associated with 
viscosity therefore is 



1 v 2 17 



Kis 30 M 2 16 



(D.41) 
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D.2 Random Statistical Fluctuations 



In the example above we have used ensemble averaging and the ergodic theorem to estimate 
average quantities in domains of size A" 1 . For any specific realization of the universe there will 
be a random statistical error in that estimate. In this section we estimate the size of that effect 
in perturbation theory. 

Stochastic Pressure 

As an illustrative example, we consider the effective pressure 

PeS=PeS + 5pi s + 5 P fi t . (D.42) 

Here, p e g denotes the renormalization of the background pressure (§4.4); bp s is the perturba- 
tion induced by long-wavelength modes 5g (§5.2); bpf^ is the perturbation induced by random 
statistical fluctuations. We will now establish the following two facts: 

1. Both 8p eS and Sp^^ are much smaller than p e g if the momentum integral defining p e g is 
UV-dominated (i.e. dominated by scales with q+ S> A). 

2. Typically, Sp^ff* > <5Pcff with equality holding at q* ~ /cnl (in Einstein-de Sitter) or q* ~ A; cq 
(in our universe). However, the fluctuations 5p^g are uncorrelated with the long- wavelength 
density fluctuations bg so their effect on the power spectrum of density fluctuations is 
subdominant for scales with k <C q*. 



Proof: 

In perturbation theory we have 



8iTGa 2 nj = y^,i<A,j - 4>,k<f>,kSij ■ (D.43) 
Consider the pressure source term, r = m, 

S = 8irGa 2 t = \(t>,i4>,i ■ (D.44) 
Recall that the spatial average over a domain of size A -1 is 

[S] A = [ W A (\x-x'\)S(x')= [ f e - i ^ + ^- x W A (\q 1 + qi \)1^3l <j> q J q ^ (D.45) 

Jx' J q 1 J q 2 " 

and its ensemble average is 

([S]a) = I J d\nq q 2 A 2 (q) = 2^Ga 2 p cS . (D.46) 
The ensemble average approximates the true average of a given realization with some statistical error 

[5] A = ([S]a> ± A[5] A . (D.47) 
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How different is [S] A from its ensemble average ([<S] A ) in each cube? To assess this we compute the variance 
<P]a[S]a> =1111 e- t(qi+q2+q3+q * yx W A (\ qi + q 2 \)W A (\q 3 + q 4 \) 

Jq 1 Jq 2 Jq 3 J q 4 



Using Wick's theorem to expand the four-point function for Gaussian fields, we find 

<[S]a[S]a)-([S]a) 2 = 2 / / W A (\q 1+ q 2 \f (gl q g2)2 ^(gi)^(g 2 ). (D.49) 



Assuming that the integrals are dominated by modes with (71,(72 3> A, we may make the following substi- 
tution, W A (\q 1 + q 2 \) w 5d(|Qi + Q 2 I/ A )- Tnis g ives 

A 3 



([5] A [5] A ) - <[S] A ) 2 « / dIn «« A *(«) ' ( D - 5 °) 

(D.51) 



and 

([S] A [S] A ) - ([S] A r _ f A 3 
<[S] A } 2 ~ 2n ql ' 

where 

/ d\nxx&.\{x) q 

f = TTT\ 2A2/ M2 ' with x = —■ (D.52) 

[J dlnxx 2 A^(a;)] 2 

Here, we have factored out the scale q+ that gives the dominant contribution in the integrals. In a matter- 
only universe q± <~ qnl, while in a universe with matter and radiation q* ~ g eq . With this normalization 
the factor / is of order unity and 

([5]a[5] a )-([5] a ) 2 A 3 

— (Pup w (D - 53) 

This implies, 

5pSff* = apeff , (D.54) 

where a is a random variable with variance 

±1 - £ • (D.55) 

Eqn. (D.54) characterizes stochastic pressure perturbations. Since we have assumed A <C <?*, this shows 

that in this case Spf^ <C p e s ■ QED ■ 



The A/q+ suppression in Eqn. (D.54) may be understood as the standard ^= suppression of 
random fluctuations, where N is the number of space domains sampled. 



Simple Estimates for Einstein-de Sitter 

We now give simple estimates of the relative sizes of the above effects in the case of an Einstein- 
de Sitter universe for which the transfer function for the Newtonian potential is trivial and 
q+ ~ (/nl- Generalizing this to the case of a universe with radiation would be straightforward. 
Specifically, we will show that the fluctuations Spfg 1 are uncorrelated with the long-wavelength 
density fluctuations bi so their effect on the power spectrum of density fluctuations is subdominant 
for scales with k <C 
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Proof: 

We capture fluctuations at a scale k by taking the smoothing scale A — > k. First, we compare the 
relative sizes of the two pressure perturbations 

5pt s ~c 2 s p5 e and Spff ~ c 2 p a . (D.56) 

We find 

M^.M (D57 ) 

|^stat|2 A 2 (fc) - ^ • ^ 

We see that the two terms are equal at k <~ ^nl, where A^nl) = A a (<7NL) = 1- For fc <C <znl we have 

Near the horizon scale Spfg*" is therefore larger than Spl s by an order of magnitude. 

However, Spfg 1 is uncorrelated with long-wavelength density perturbations St so its effect on the 
density power spectrum is in fact suppressed. To see this consider the evolution equation for density 
perturbations 

S e +n5 £ + ln 2 5 e = -^M (D .59) 
2 p 

= k 2 c 2 s [S t + a] . (D.60) 
Using S for the solution in the c s = limit, we estimate the pressure-induced corrections 

k 2 k 2 

Si « St + + ^c 2 s a . (D.61) 

The first correction term (arising from tidal pressure) is correlated with St, while the second (arising from 
stochastic pressure) is not. Hence, the corrections to the power spectrum are 

^ = AI+^A?+(^Va^ (D.62) 



and 



Using 



we find 



AP| fk 2 ,\-A 2 



APf at \H 2 s 1 A 2 





k 2 






U 2 










c 

H 2 




\9nl 


AP| 




<?NL 



(D.63) 



c 2 s ~ |^A^( 9NL ) and £cj - ( " ) A^ NL ) = ('^ )' A^-) (D.(S4) 



AP| tat fe 



(D.65) 



We see that the corrections to the power spectrum are equal at the non-linear scale, but dominated by 
tidal pressure on larger scales, k~ x ^> g^ L . QED ■ 
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